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Distributions  of  Magnetic  Field 
Variations,  Differences  and  Residuals 

J.  Watermann,  J.  Lam 


Executive  Summary: 

The  search  for  submarines  by  means  of  Magnetic  Anomaly  Detection  (MAD)  is 
effected  by  a  low-flying  aircraft  with  a  magnetometer  installed  onboard,  or,  for 
better  sensitivity,  mounted  on  an  external  boom.  The  magnetic  field  produced 
by  a  submarine  appefirs  as  a  local  quasistatic  distortion  (’’anomaly”)  of  the 
geomagnetic  field.  The  technique  requires  a  remote  reference  magnetometer  at 
a  fixed  place  in  order  to  compensate  for  the  omni-present  temporal  variations 
of  the  geomagnetic  field.  This  place  may  be  many  kilometres  away  from  the 
position  of  the  aircraft,  in  which  case  the  technique  requires  a  high  spatial 
uniformity  of  the  ambient  magnetic  field  variations. 

While  application  of  the  MAD  technique  may  be  possible  in  many  areas,  there 
are  some  areas  where  this  method  must  fail.  These  are  either  areas  with  poorly 
mapped  magnetostatic  anomalies  of  small  spatial  scale,  or  areas  with  high 
levels  of  non-imiform,  man-made  magnetic  noise.  The  latter  typically  include 
coastal  zones  with  industrial  and  service  centers  in  the  vicinity.  We  find  among 
them  the  southern  Ligurian  coastal  zone  which  is  magnetically  highly  disturbed 
by  a  busy  electrified  railway  line.  Electric  train  supply  currents  leak  several 
kilometres  out  into  the  sea  and  produce  non-uniform  magnetic  field  distortions 
(anomalies)  of  large  amplitude. 

We  propose  to  employ  in  such  an  area  the  Moving  Target  Detection  (MTD) 
technique.  Compared  to  MAD,  the  roles  of  sensor  and  submarine  are  inverted. 
Several  magnetometers  reside  in  fixed  positions  on  the  sea  bottom  at  suitable 
array  nodes  and  detect  localized  temporal  magnetic  field  distortions  (such  as 
those  created  by  moving  submarines)  through  comparison  between  simultane¬ 
ous  measurements  from  different  sensors.  Depending  on  the  spatial  uniformity 
of  the  ambient  magnetic  field,  the  sensor  spacing  may  be  large  or  must  be  small. 

We  found  that  in  an  area  some  50  km  from  the  coast,  a  magnetometer  located 
several  kilometres  away  from  the  zone  kept  under  surveillance,  can  provide  for 
a  good  magnetic  field  reference  and  is  not  influenced  by  the  magnetic  field 
of  a  submarine  navigating  in  the  surveillance  area.  In  the  southern  Ligurian 
coastal  zone,  reference  magnetometers  must  be  very  close  to  the  area  under 
surveillance  (one  kilometre  or  less  apart)  in  order  to  provide  for  satisfactory 
magnetic  field  compensation.  In  such  a  case,  the  reference  sensor  may  also  be 
influenced  by  the  magnetic  field  of  the  submarine,  and  surveillance  by  means 
of  magnetic  field  observations  is  more  complicated.  It  is  nevertheless  possible, 
even  with  limited  computing  power. 
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We  found  that  the  vertical  component  of  the  ambient  magnetic  field  is  difficult 
to  compensate  for;  it  contains  large  signals  not  linearly  correlated  between  ad¬ 
jacent  sites  and  can  not  serve  as  a  suitable  detector.  The  smallest  uncorrelated 
contributions  to  the  magnetic  field  variations  were  found  when  observations  of 
the  coast^perpendicular  magnetic  field  component  from  sites  at  equal  distance 
from  the  coast  were  compared.  Surveillance  arrays  should  therefore  be  sited  to 
form  pairs  with  magnetometers  at  equal  distance  from  the  coast. 

For  the  type  of  magnetometers  used  in  our  experiment  and  for  arbitrarily  se¬ 
lected  magnetic  field  perturbation  thresholds,  we  give  the  rate  at  which  the 
specified  threshold  is  likely  to  be  exceeded,  due  to  signals  from  natural  and 
man-made  sources  other  than  surface  vessels  or  submarines.  This  rate  can  be 
considered  the  false  alarm  rate  of  the  magnetic  surveillance  system.  To  reduce 
the  false  alarm  rate,  the  surveillance  system  operator  must  increase  the  detec¬ 
tion  threshold  or  reduce  the  magnetometer  spacing.  The  report  concludes  with 
an  example  in  which  we  specify  a  threshold  of  0.4  nT  and  1  km  sensor  spac¬ 
ing.  We  show  that  several  false  alarms  per  hour  would  be  issued  if  the  vertical 
magnetic  field  component  were  considered,  but  practically  no  false  alarm  if  the 
coast-perpendicular  component  of  a  pair  of  magnetometers  at  equal  distance 
from  the  coast  were  considered. 

One  of  the  critical  issues  in  designing  a  magnetic  surveillance  system  is  the  se¬ 
lection  of  the  magnetometer  type.  Our  observations  in  the  coastal  zone  suggest 
that  the  spatial  non-uniformity  of  the  magnetic  field  requires  sensors  with  small 
spacing  (not  more  than  one  kilometre)  for  satisfactory  ambient  magnetic  field 
compensation.  This  implies  in  many  cases  a  large  number  of  sensors.  Sophis¬ 
ticated  instruments  with  extremely  low  noise  and  high  stability  characteristics 
are  less  useful.  Compensation  works  best  if  all  three  vector  components  of  the 
magnetic  field  are  used.  This  requires  the  use  of  tri-axial  sensors  rather  than 
scalar  magnetometers. 

We  need  to  emphasize  that  these  results  are  area-specific.  Once  an  area  has 
been  selected  for  magnetic  surveillance,  magnetic  field  measurements  using  at 
least  two  spaced  sensors  must  be  made  prior  to  deploying  the  surveillance  sys¬ 
tem,  in  order  to  determine  which  arrangement  of  magnetometers  and  which 
magnetic  field  component  is  best  suited  for  surveillance  purposes.  Such  a 
preparatory  survey  would  typically  require  a  few  days  of  measurements  (at 
least  one  or  two  full  days),  preferably  continuously.  Due  to  the  physical  prop¬ 
erties  of  the  ambient  magnetic  field  close  to  the  coast,  the  vertical  component 
will  probably  never  be  very  useful.  We  also  expect  that  the  strike  direction  of 
the  coast  will  usually  favour  an  arrangement  which  builds  on  the  correlation 
between  sensors  at  the  same  distance  from  the  coast. 

This  report  demonstrates  that  the  efficiency  of  a  magnetic  surveillance  system 
using  an  array  of  fixed  sensors  (which  we  call  the  MTD  technique)  depends 
largely  on  the  arrangement  of  the  sensors  relative  to  the  coast  and  on  the 
magnetic  field  vector  component  selected  for  detecting  a  passing  subm2u:ine. 
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Distributions  of  Magnetic  Field 
Variations,  Differences  and  Residuals 

J.  Watermann,  J.  Lam 


Abstract: 


Temporal  and  spaticil  variations  of  the  geomagnetic  field  were  recorded  in  differ¬ 
ent  geographic  areas,  using  arrays  of  seven  tri-axial  magnetometers  on  the  sea 
bottom  at  50-150  m  depth  for  periods  of  up  to  16  days.  The  different  geographic 
areas  are  characterized  by  different  levels  of  the  mean  ambient  magnetic  noise. 
We  discuss  data  from  a  magnetically  quiet  area  and  then  focus  on  a  coastal 
zone  highly  disturbed  by  anthropogenic  magnetic  fields,  in  particular  by  noise 
from  an  electrified  railway  line.  In  the  latter  area,  the  spatial  uniformity  of 
the  ambient  magnetic  field  is  rather  poor  and  the  correlation  between  adjacent 
sites  lower  than  in  the  former  one. 


Simultaneous  magnetic  field  observations  from  pairs  of  neighbouring  magne¬ 
tometers  are  analysed  and  compared.  We  compute  first  vector  differences  of 
the  magnetic  variations  measured  at  neighbouring  sites,  and  subsequently  vec¬ 
tor  residuals.  The  residuals  are  those  contributions  to  the  magnetic  field  vari¬ 
ations  which  are  not  correlated  between  adjacent  sites  and  which  can  not  be 
represented  through  a  linear  tri variate  model. 

In  the  magnetically  quiet  area,  we  find  that  the  magnetic  field  residuals  are  of 
the  order  of  the  system  noise  for  magnetometers  with  about  1  km  spacing  and 
slightly  higher  for  those  with  12  km  spacing.  The  area  lends  support  to  the 
application  of  a  Remote  Reference  Technique,  i.e.  a  technique  in  which  the  am¬ 
bient  magnetic  field  is  compensated  for  by  using  measurements  from  a  remotely 
operated  reference  sensor.  "Remote”  means  that  the  reference  sensor  is  suffi¬ 
ciently  distant  not  to  be  influenced  by  the  magnetic  field  of  a  ship  navigating 
in  the  vicinity  of  the  surveillance  magnetometers. 


Although  the  magnetic  field  variations  in  the  highly  disturbed  coastal  zone 
are  several  times  greater  than  typical  magnetic  field  variations  of  natural  ori¬ 
gin,  and  also  spatially  non-uniform,  we  find  good  linear  correlation  between 
adjacent  sites.  This  is  demonstrated  by  the  often  rather  small  magnetic  field 
residuals.  However,  in  order  to  achieve  small  residuals,  the  spacing  between  the 
magnetometers  must  be  an  order  of  magnitude  smaller  than  in  the  quiet  zone. 
This  may  require  the  use  of  a  Local  Reference  Technique  for  compensation  of 
the  ambient  magnetic  field.  In  a  Local  Reference  Technique,  the  reference  sen¬ 
sor  must  be  so  close  to  the  surveillance  sensors  that  it,  too,  is  influenced  by  the 
magnetic  field  of  the  submarine  navigating  in  the  vicinity  of  the  surveillance 
sensors. 
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In  the  coastal  zone^  we  find  that  the  residuals  of  the  vertical  component  of  the 
magnetic  field  variations  remain  rather  large.  The  coast-perpendicular  com¬ 
ponent  of  magnetometers  at  equal  distance  from  the  coast  yields  the  smallest 
residuals,  i.e.  they  yield  the  best  compensation  of  magnetic  field  variations 
when  using  a  reference  technique.  The  intensity  distribution  of  magnetic  field 
variations  observed  at  the  various  sites  strongly  suggests  that  a  significant  part 
of  the  electric  leakage  from  the  railway  system  flows  out  into  the  sea  up  to 
several  kilometres  from  shore. 

In  summary,  we  find  that  magnetic  field  variations  in  the  southern  Ligurian 
coastal  zone  are  dominated  by  magnetic  railway  noise  which  results  in  an 
anisotropic  distribution  of  the  spatial  correlation  between  spaced  sites. 
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Introduction 


The  steel  mass  of  a  ship  provides  for  a  magnetic  field  which  is  superimposed  on 
the  existing  ambient  magnetic  field  and  modifies  it  locally.  This  modification  is  de¬ 
tectable  with  commercially  available  magnetometers  at  ranges  of  hundreds  of  metres 
and  sometimes  up  to  a  few  kilometres.  A  ship’s  magnetic  moment  consists  of  two 
components,  a  permanent  one  and  an  induced  one.  With  respect  to  their  magneti¬ 
zation  and,  consequently,  their  magnetic  signature,  there  is  no  difference  in  principle 
between  surface  vessels  and  submarines. 

The  permanent  magnetic  moment  is  impressed  on  the  steel  in  the  production  process, 
at  the  time  when  its  temperature  falls  below  the  Curie  point.  "Permanent”  in 
the  sense  used  here  does  not  imply  that  it  absolutely  never  changes.  However, 
permanently  magnetized  material  needs  to  be  exposed  to  a  ojnstant  magnetic  field 
over  a  long  time  (several  months,  and  sometimes  even  years)  in  order  to  acquire, 
at  least  partially,  a  magnetic  moment  having  the  orientation  of  the  applied  ambient 
field  and  thus  changing  the  material’s  effective  moment.  The  same  holds  for  the 
decay  of  the  permanent  magnetic  moment. 

The  induced  magnetic  moment  varies  instantaneously  with  the  intensity  and  orien¬ 
tation  of  the  ambient  field  to  which  the  material  is  exposed.  In  the  open  sea,  the 
dominant  ambient  magnetic  field  is  the  geomagnetic  field,  99.9%  of  which  stem  firom 
quasi-permanent  internal  sources  (their  variations  are  usually  called  "secular”  be¬ 
cause  they  become  significant  only  over  a  period  of  many  years).  This  implies  that 
the  induced  magnetic  moment  changes  with  the  position  of  the  vessel  with  respect 
to  its  "geomagnetic  coordinates”. 

"Geomagnetic  coordinates”  is  used  here  as  a  loosely  defined,  generic  name.  Many 
different  geomagnetic  coordinate  systems  with  different  properties  exist.  For  the 
purpose  of  getting  a  sufficiently  accurate  quantitative  representation  of  the  geomag¬ 
netic  field  at  any  given  geographic  position,  altitude  above  sea  level,  and  geomagnetic 
epoch,  the  International  Geomagnetic  Reference  Field  (IGRF)  [1]  is  a  good  source 
to  start  with.  It  represents  the  internal  global  component  of  the  geomagnetic  field 
and  is  updated  every  five  years.  The  presently  valid  version,  the  IGRF-95,  covers 
the  years  1945  through  2000.  The  IGRF  uses  a  spherical  harmonic  expansion  up 
to  degree  and  order  10  and  can  therefore  not  very  well  represent  regional  and  lo¬ 
cal  magnetic  anomalies.  Regional  geomagnetic  field  models  have  been  compiled  for 
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various  parts  of  the  earth,  including  the  Central  Mediterranean  area  (i.e.  Italy  and 
surrounding  seas),  for  which  an  expansion  in  spatial  spherical  harmonics  of  degree 
four  and  order  two  and  in  temporal  polynomials  of  degree  three  exists  [2] .  Regional 
geomagnetic  field  models  for  Germany  axe  published  by  the  ministry  of  transport  [3]. 
The  local  magnetic  field  can  usually  be  taken  fi:om  anomaly  maps  where  available, 
with  reasonably  accuracy  and  resolution  in  the  vicinity  of  land  masses  where  the 
density  of  measiurements  is  often  high,  less  accmate  in  the  open  sea  where  measure¬ 
ments  are  sparse.  For  Italy  and  the  Ligmrian,  Tyrrhenian,  Ionic  and  Adriatic  Seas, 
magnetic  anomaly  maps  are  periodically  published,  e.g.  [4]. 

If  the  internal  sources  were  the  only  ones  which  determine  the  intensity  and  ori¬ 
entation  of  the  geomagnetic  field,  the  problem  of  detecting  a  steel-hulled  vessel  by 
observing  its  magnetic  signature  were  reduced  to  determining  as  accurately  as  pos¬ 
sible  the  ambient  magnetic  field  at  a  given  location  in  absence  of  ships,  and  then  use 
this  information  as  the  standard  against  which  subsequent  magnetic  field  measure¬ 
ments  (at  the  same  location)  are  to  be  compared.  Any  deviation  from  the  standard 
would  then  be  explained  by  the  presence  of  a  magnetic  source  subsequently  intro¬ 
duced  into  the  area,  such  as  a  ship  bearing  a  magnetic  moment.  In  reality,  the 
ambient  magnetic  field  in  a  marine  environment  is  composed  of  contributions  firom 
various  types  of  magnetic  field  sources:  (a)  quasi-permanent  soiurces  of  internal  and 
anthropogenic  origin;  (b)  oceanic  sources,  moderately  variable  in  space  and  time; 
(c)  extraterrestrial  sources,  highly  variable  in  time,  much  less  variable  in  space;  (d) 
anthropogenic  sources,  highly  variable  in  space  and  time. 

Two  principally  diflFerent  magnetic  detection  techniques  are  nowadays  in  use.  One 
is  known  by  the  acronym  MAD  (Magnetic  Anomaly  Detection).  In  this  technique, 
large  areas  of  the  sea  are  magnetically  surveyed  by  planes.  If  a  submarine  is  present, 
its  magnetic  field  is  observed  as  a  local  deviation  of  the  tqtaJ  magnetic  field  from  its 
undisturbed  level,  i.e.  the  submarine  is  considered  a  local  magnetic  anomaly.  If  the 
existing  magnetic  anomalies  such  as  ship  wrecks,  pipelines,  power  cables,  or  geologi¬ 
cally  complex  sea  bed  structures,  are  known,  or  if  the  seh  is  sufficiently  deep  so  that 
their  effect  can  not  be  observed  at  the  surface,  any  previously  unknown  anomaly  can 
be  identified  and  possibly  related  to  a  submarine  if  no  surface  vessel  is  observed.  Be¬ 
cause  the  planes  fly  many  times  faster  than  submarines  move,  this  technique  means 
essentially  searching  with  a  moving  magnetometer  for  a  static  magnetic  anomaly. 

The  other  technique  which  is  explored  at  SACLANTCEN,  inverts  the  roles  of  sen¬ 
sor  and  target.  A  magnetometer  is  installed  on  the  sea  bottom,  fixed  in  location 
and  attitude,  and  records  temporal  variations  of  the  ambient  magnetic  field.  Those 
magnetic  anomalies  which  remain  fixed  in  place  are  not  noticed  by  the  magnetome¬ 
ter.  Only  a  moving  target,  e.g.  a  moving  submarine,  which  represents  a  moving 
magnetic  moment,  generates  a  time  dependent  magnetic  field  and  leaves  a  signal  at 
the  magnetometer.  We  will  use  the  acronym  MTD  (Moving  Target  Detection)  for 
this  technique.  Ambient  magnetic  field  sources  of  type  (a)  from  the  list  above  play 
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no  role  in  the  MTD  technique  but  they  may  have  a  significant  impact  on  the  MAD 
technique. 

So  far,  the  MAD  technique  employs  in  practice  sensors  which  measiure  only  the 
magnetic  field  intensity  but  not  its  orientation.  This  restriction  is  imposed  by  a 
sensor  attitude  control  problem.  A  rotation  of  a  vector  magnetometer  by  just  four 
arcseconds  can  have  the  same  effect  on  the  measurements  as  a  change  of  the  ambient 
magnetic  field  by  1  nT  (nanoTesla).  For  comparison,  a  typical  1500-ton  submarine 
not  treated  for  magnetic  compensation  would  generate  a  magnetic  field  of  1  nT  at 
some  300  m  range.  A  sensor  fixed  on  the  sea  bottom  in  a  stable  position  can  take 
advantage  of  measuring  the  magnetic  field  variation  in  all  three  vector  components. 
Here,  only  sensor  oscillations  forced  by  the  moving  sea  water  can  pose  a  problem  if 
insufficient  care  is  taken  [5]. 

The  positioning  problem  of  the  MAD  technique  is  less  grave  than  the  attitude  prob¬ 
lem,  but  it  exists.  A  plane  flying  over  the  Tyrrhenian  Sea  would  measture  a  change 
of  1  nT  by  changing  its  altitude  by  50  m,  or  its  position  by  500  m  in  north-south 
direction.  South  of  Montecristo,  in  an  area  characterized  by  an  intense  natural  mag¬ 
netic  field  anomaly,  a  horizontal  displacement  by  just  30  m  is  sufiicient  to  chaJige 
the  observed  magnetic  field  by  1  nT.  Obviously,  this  problem  does  not  exist  for  a 
magnetometer  residing  at  a  fixed  place  on  the  sea  bottom. 

Potential  problems  common  to  both  techniques  are  listed  in  items  (b)-(d)  above. 
For  the  areas  investigated  in  the  SACLANTCEN  programme  on  environmental  geo¬ 
magnetic  measurements,  it  could  be  demonstrated  in  [5]  that  the  magnetic  dynamo 
field  from  sea  water  motion  is  detectable  but  of  such  small  intensity  that  it  is  neg¬ 
ligible  in  comparison  with  the  other  sources  and  thus  plays  no  role  for  the  MAD 
and  MTD  techniques.  Only  items  (c)  and  (d),  extraterrestrial  and  variable  anthro¬ 
pogenic  sources,  remain  to  cause  problems.  They  are  the  subject  of  our  report. 
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Data  Coverage 


We  discuss  in  this  report  magnetic  field  measurements  made  on  the  sea  bottom  along 
the  southern  Ligurian  coast  (Cruise  1996)  and  near  Isola  di  Montecristo  (Cruise 
1997).  Area  maps  are  given  below  (Maps  1  and  2).  Tables  1  and  2  list  the  WGS-84 
coordinates,  bottom  depths,  and  data  recording  intervals  of  those  magnetometer 
modules  which  were  operated  and  recovered  successfully.  Our  report  demonstrates 
that  the  near-coast  array  recorded  highly  disturbed  spatially  varying  signals  re¬ 
sulting  firom  the  proximity  of  an  active  railway  line  (Cruise  1996).  In  contrast, 
measurements  made  some  50  kilometres  from  the  coast  in  an  electromagnetically 
very  quiet  environment  display  a  high  degree  of  spatial  uniformity  (Cruise  1997). 


Table  1  1996  Measurements  (see  also  Map  1) 


Module 

WGS>84  coordinates 

Bottom  depth 

Data  interval  (UTC) 

M03 

44°  15.76'  N 

009°15.47'  E 

116  m 

10/08,  14:00-10/14,  11:45 

M04 

44°05.79'  N 

009°35.73'  E 

130  m 

09/28,  20:00-10/14, 11:45 

M06 

44°08.20'  N 

009°30.79'  E 

127  m 

09/28,  20:00-10/14,  11:45 

M12 

44°  10.09'  N 

009°27.15'  E 

132  m 

09/28,  20:00-10/14, 11:45 

M17 

44°  12.70'  N 

009°29.72'  E 

58  m 

10/04,  10:00-10/14, 11:45 

M18 

44°11.98'  N 

009°29.12'  E 

83  m 

09/28,  20:00-10/14, 11:45 

M19 

44°11.11'  N 

009°28.26'  E 

108  m 

09/28,  20:00-10/14, 11:45 

We  used  in  both  sea  trials  the  same  type  of  autonomous  Ocean  Bottom  Magneto¬ 
meters  (OEMs).  They  are  complex  modules  which  integrate  a  tri-axial  linear-core 
AC-coupled  fluxgate  magnetometer,  two  ±25°  inclinometers  aligned  with  the  x-  and 
y-axes  of  the  magnetic  field  sensor,  and  a  2-D  fiuxgate  compass  mounted  parallel  to 
the  x-y-plane  which  serves  to  identify  the  magnetic  bearing  of  the  sensor’s  a:-axis 
once  the  module  has  settled  on  the  sea  bottom.  A  crystal  oscillator  with  ±1  ppm 
long-term  stability  controls  the  digital  data  acquisition  system.  The  manufacturer’s 
design  of  the  fiuxgate  feedback  circuit  and  our  setting  of  the  operating  parameters 
result  in  a  sampling  rate  of  2  Hz  and  a  usable  frequency  bandwidth  of  1-500  mHz. 
The  data  are  recorded  in  16-bit  resolution,  with  the  least  significant  bit  equaling 
3.05  pT  (picoTesla),  and  stored  in  non-volatile  solid  state  memory  (fiash  memory 
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Table  2  1997  Measurements  (see  also  Map  2) 


Module 

WGS-84  coordinates 

Bottom  depth 

Data  interval  (UTC) 

M03 

42°21.88'  N 

010°18.70'  E 

159  m 

07/18,  09:30-08/02,  01:00 

M04 

42°21.82'  N 

010°18.70'  E 

159  m 

07/18,  09:00-07/26, 11:45 

M06 

42°21.66'  N 

010n8.70'  E 

157  m 

acuisition  failed 

M12 

42°21.55'  N 

010°18.71'  E 

144  m 

very  noisy  30%  of  the  time 

M17 

42°23.95'  N 

010°10.50'  E 

134  m 

07/19,  14:30-07/26, 13:15 

M18 

42°24.15'  N 

010°  10.50'  E 

134  m 

07/19,  15:00-07/26, 13:30 

M19 

42°24.55'  N 

010°  10.50'  E 

132  m 

07/19,  15:30-07/26, 13:45 

caxds).  Upon  retrieval  of  the  modules,  the  data  are  transferred  from  the  flash  cards 
to  an  IBM  compatible  PC  and  subsequently  written  to  CD-ROM. 
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The  Submarine  Magnetic  Field 


The  magnetization  of  a  submarine  can  be  represented  by  a  series  of  magnetic  mul¬ 
tipoles  with  common  origin.  According  to  widely  accepted  theory  (V  •  B  =  0),  no 
magnetic  monopoles  exist,  therefore  the  dipole  is  the  lowest  order  contribution  to 
the  multipole  expansion  series.  The  next  higher  multipole,  the  quaxirupole,  can  be 
eliminated  from  the  sphericcil  harmonics  expansion  if  the  z-axis  of  the  reference  co¬ 
ordinate  system  is  aligned  with  the  dipole  moment  and  the  location  of  the  magnetic 
origin  properly  selected  [6].  The  submarine’s  center  of  gravity,  for  instance,  would 
not  normally  be  a  good  choice  for  its  magnetic  origin.  The  submarine’s  magnetic 
field,  can  thus  be  written  in  the  form 


g(sm) 


Air 


3(M-f)f-M  ^/IM 


(1) 


M  denotes  the  dipole  moment,  r  the  distance  vector  from  the  submarine’s  magnetic 
center  to  the  magnetometer  (f  the  corresponding  unit  length  vector),  and  0(f) 
means  that  the  rest  of  the  series  is  asymptotically  limited  by  the  fimction  /  times  a 
constant  factor.  Obviously,  the  effect  of  the  octopole  and  higher  multipoles  decreases 
rapidly  with  increasing  distance  between  submarine  aiid  magnetometer.  A  rule  of 
thumb  which  was  derived  from  measurements  on  cars  and  ships  suggests  that  for  a 
distance  in  excess  of  about  1.5  times  the  length  of  the  object,  only  the  dipole  moment 
remains  to  be  of  significance.  An  accurate  quantitative  statement,  however,  depends 
on  the  magnetization  of  the  individual  object  and  would  need  to  be  confirmed  by 
specific  measurements.  Such  measmements  have  not  been  made  at  SACLANTCEN. 
Magnetic  signatures  from  USN  submarines,  collected  by  NRaD  and  analyzed  by 
Hughes  Aircraft  [7],  indicate  that  the  higher  multipoles  may  have  some  bearing  on 
the  composition  of  those  signals  which  were  observed  at  close  distance  (some  100 
feet),  but  quantitative  results  were  not  given.  It  is  also  not  the  purpose  of  this  report 
to  comment  quantitatively  on  the  relative  importance  of  the  magnetic  octopole  and 
higher  multipoles. 

Let  us  assume  that  the  magnetometer  is  always  sufficiently  far  from  the  submarine 
so  that  it  observes  only  its  dipole  field  while  the  field  from  the  higher  multipoles 
is  absorbed  in  the  system  noise.  We  can  then  represent  each  of  the  three  vector 
components  of  the  magnetic  perturbation  generated  by  a  submarine  sailing  along 
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a  straight  line  with  constant  speed,  through  a  linear  combination  of  three  simple 
functions  named  ’’Anderson  functions”: 


Ao(r)  =  (t2  +  1)-5/2 

Ai(r)  =  r(T2 +  l)-5/2 
A2{t)  =  r2(T2  +  l)-5/2 


The  dimensionless  parameter 

r  =  — 
ro 


(2) 

(3) 


relates  submarine  speed,  vq,  and  distance  of  closest  approach,  tq,  to  the  elapsed 
time,  t,  centered  on  the  time  of  closest  approach  [8]. 

Fig.  1  is  a  graphical  representation  of  the  Anderson  functions,  using  the  parameters, 
submarine  speed,  6  knts,  and  distance  at  closest  approach,  300  m.  Theoretically, 
the  magnetic  field  variation  seen  by  each  of  the  three  axes  of  a  vector  magnetometer 
consists  of  a  linear  combination  of  these  three  functions,  with  coefficients  depending 
on  the  orientation  of  the  dipole  moment  relative  to  the  major  axis  of  the  submarine 
and  the  submarine’s  coxirse.  When  looking  for  the  magnetic  signature  of  a  submarine 
one  needs  to  look  for  magnetic  field  transients  composed  of  the  Anderson  functions. 

We  have  no  reason  to  assume  that  the  dipole  moment  of  a  ship  is  oriented  parallel 
to  its  major  axis  (i.e.  along  a  horizontal  line  from  bow  to  stern).  It  can,  in  fact, 
have  any  arbitrary  direction,  depending  on  the  history  of  the  vessel  (which  affects  its 
permanent  magnetization)  and  its  actual  location  and  orientation  (which  determines 
its  induced  magnetization. 

An  example,  based  on  measmrements  made  in  1997,  close  to  Isola  di  Montecristo,  is 
shown  in  Fig.  2  (see  Map  2  for  reference).  The  SACLANTCEN  workboat  Manning 
sailed  with  a  constant  speed  of  10.3  knts  on  a  straight  west-east  course  right  over 
OEM  04,  which  resided  1.5  km  north  of  the  island  at  159  m  depth.  The  boat  crossed 
it  at  14:38  UTC.  A  second  sensor,  OEM  03,  was  located  110  m  north  of  OEM  04,  also 
at  159  m  depth.  Only  the  vector  residuals  between  OEMs  03  and  04,  and  between 
those  two  and  a  third  sensor,  OEM  19,  located  12  km  northwest  of  the  former  two, 
eire  displayed  while  the  imcompensated  magnetic  field  variations,  which  are  an  order 
of  magnitude  larger,  axe  suppressed  here.  The  residuals  are  basically  differences 
between  magnetic  field  vectors  from  two  sensors,  computed  after  optimal  alignment 
and  linear  scaling  of  the  vector  components.  Thus,  the  residual  is  what  remains 
after  optimal  linear  magnetic  field  compensation  with  the  help  of  a  reference  sensor. 
The  term  ’’optimal  alignment  and  linear  scaling”  is  explained  more  precisely  in  a 
subsequent  section.  The  residuals  to  OEM  19  (blue  and  green  lines)  suggest  that 
Manning’s  magnetic  dipole  moment  was  confined  to  an  east-vertical  plane  and 
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pointed  almost  vertically  downward,  i.e.  it  was  not  aligned  with  the  geomagnetic 
field.  We  conclude  that  Manning’s  magnetic  dipole  moment  is  dominated  by  its 
permanent  component  while  the  induced  one  is  negligible. 

A  rule  of  thumb,  derived  from  measurements  on  various  vessels  and  submarines 
reported  in  the  classified  literature,  suggests  that  the  magnetic  moment  of  an  un¬ 
treated  ship  corresponds  to  about  400  A-m^  times  the  steel  mass  (in  tons)  of  the  ship. 
According  to  this  rule,  the  magnetic  dipole  field  of  Manning,  as  it  is  seen  by  OBM 
04,  should  amount  to  about  1.6  nT,  50%  more  than  the  1.1  nT  actually  observed. 
At  the  location  of  OBM  03  we  would  have  expected  about  0.9  nT  as  opposed  to  the 
0.7  nT  actually  observed. 

The  same  rule  of  thumb  provides  for  an  estimate  of  the  detection  range  for  a  sub¬ 
marine  by  means  of  magnetic  field  sensors.  If  the  effective  ambient  noise  plus  the 
instrument  noise  remain  below  0.25  nT  in  the  frequency  band  of  interest,  and  the 
detection  threshold  is  set  to  0.25  nT,  a  1500-ton  submarine  might  be  detected  if 
it  passes  the  sensor  within  a  range  of  some  600  m.  The  rule  of  thumb  can  also  be 
used  to  estimate  the  time  interval  a  submarine  navigates  within  the  detection  range. 
If  our  sample  submarine  sails  with  100  m  clearance  at  6  knts  speed  directly  over 
the  magnetometer,  its  magnetic  field  would  exceed  0.25  nT  for  a  period  of  about  6 
minutes. 
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Magnetic  Field  Reference  Techniques 


A  permanently  present  background  of  natural  magnetic  field  fiuctuations  covering  a 
wide  firequency  band,  including  the  firequency  band  of  interest  to  ASW,  provides  for 
a  level  of  magnetic  ambient  noise  the  intensity  of  which  is  usually  much  higher  than 
the  sensitivity  of  the  magnetometers.  Fig.  3,  compiled  fi:om  data  found  in  [9],  [10], 
shows  a  distribution  of  the  mean  peaJc-to-peak  amplitudes  of  geomagnetic  pcl-pc5 
fluctuations,  observed  over  many  years  in  Central  Europe  (more  specifically,  in  the 
vicinity  of  Gottingen,  F.R.G.).  "Mean  amplitude”  means  that  an  average  was  taken 
over  all  those  time  intervals  in  which  fluctuations  within  a  given,  narrow  frequency 
band  could  be  identified.  Consequently,  geomagnetic  ULF  fluctuations  possess  at 
times  amplitudes  exceeding  the  mean  values  displayed  here.  The  pc  type  of  ULF 
fluctuations  is  just  one  specific  type  (and  usually  not  the  most  intense  one)  of  a  wide 
range  of  different  types  of  geomagnetic  field  fluctuations,  and  is  observed  practically 
every  day,  and  often  for  many  hours. 

Fig.  1  demonstrates  that,  for  the  submarine  parameters  assumed  here,  the  magnetic 
field  pertmrbation  would  have  a  typical  dmration  of  300-400  seconds.  From  Fig.  3 
we  learn  that  geomagnetic  pc  fluctuations  in  the  corresponding  period  range  reach 
an  average  amplitude  of  6  nT  p-p.  Such  an  amplitude  would  be  equivalent  to  a 
magnetic  signal  firom  our  sample  submarine  navigating  at  less  than  200  m  range.  In 
order  to  increase  the  detection  range,  a  technique  has  to  be  developed  which  aids  in 
suppressing  the  effect  of  the  geomagnetic  fiuctuations. 

Building  on  the  fact  that  the  most  important  contribution  to  the  natural  geomag¬ 
netic  fluctuations  stems  firom  sources  in  the  ionosphere  and  magnetosphere,  more 
than  100  km  above  the  ground,  and  the  second  important  contribution  from  non- 
imiform  conductivity  distributions  found  at  geographic  discontinuities  such  as  coast 
lines,  we  can  assume  that  in  most  areas  the  spatial  uniformity  of  the  natural  geo¬ 
magnetic  fluctuations  is  sufficiently  high  over  a  distance  of  at  least  a  few  kilometres 
to  permit  simple  linear  magnetic  fleld  compensation.  A  more  detailed  discussion 
follows  further  below,  along  with  the  presentation  of  magnetic  field  data  from  the 
1997  cruise. 

The  spatial  uniformity  of  the  geomagnetic  field  fluctuations  is  employed  in  compiling 
the  example  shown  in  Fig.  4.  Fig.  4a  displays  vector  time  series  firom  two  different 
events  recorded  by  OBM  04  in  September  1995  in  the  vicinity  of  the  Formiche  di 
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Grosseto  group  of  islands  (see  Map  3  for  reference).  In  each  panel,  a  section  is 
highlighted  which  could,  according  to  Fig.  1,  possibly  be  associated  with  a  passing 
vessel.  Fig.  4b  displays  a  stack  of  simultaneous  measiurements  from  three  magne¬ 
tometers  (eastward  components  only),  spanning  a  distance  of  eight  kilometres  (c.f. 
Map  3).  The  upper  panel  of  Fig.  4b  reveals  that  the  suspicious  signal  observed  on 
September  16  around  18:41  UTC  appeared  at  all  three  sites  with  almost  identical 
amplitude  and  shape.  The  mean  difference  field  (the  difference  between  OBM  04 
and  a  weighted  mean  of  OEMs  18  and  19)  is  not  pjirticularly  enhanced.  A  ship 
signature  can  therefore  be  excluded.  At  8  km  distance,  even  a  ship  with  100,000 
tons  of  steel  mass  would  generate  a  magnetic  field  of  a  couple  of  pT  (picoTesla) 
only,  much  less  than  the  system  noise  of  our  magnetometers.  The  lower  panel,  how¬ 
ever,  demonstrates  that  the  suspicious  signal  observed  oh  September  20  at  08:46/:47 
UTC  was  seen  only  by  OBM  04  and  not  by  the  OEMs  18  and  19.  It  must  have  been 
generated  by  a  local  source,  for  which  we  suggest  a  passing  vessel. 

Fig.  4  demonstrates  that  fluctuations  of  the  ambient  magnetic  field  can  be  compen¬ 
sated  for  by  using  a  reference  magnetometer  as  long  as  its  distance  fi-om  the  surveil¬ 
lance  magnetometer  (termed  "test  magnetometer”  in  the  following)  is  significantly 
smaller  than  the  spatial  scale  length  of  the  ambient  magnetic  field  fluctuations.  This 
may  in  some  cases  permit  to  choose  a  distance  which  exceeds  the  range  within  which 
the  ship  generates  a  measurable  signal.  In  such  a  case,  when  the  test  magnetometer 
records  a  magnetic  signal  from  a  ship,  the  reference  magnetometer  is  free  of  such 
a  signal  and  thus  provides  for  a  proper  magnetic  field  baseline.  We  term  this  the 
Remote  Reference  Technique.  Fig.  2  may  serve  as  an  example  for  this  technique 
(blue  and  green  lines).  OEMs  03  and  04  record  the  magnetic  signature  of  the  work- 
boat  Manning  while  OBM  19,  12  km  away,  is  guaranteed  free  from  any  magnetic 
influence  by  Manning.  The  residual  shows  Manning’s  true  magnetic  signature.  It 
appears  of  larger  amplitude  at  OBM  04,  and  of  smaller  but  still  detectable  amplitude 
at  OBM  03. 

In  other  cases,  when  the  ambient  magnetic  field  is  l^s  uniform,  a  satisfactory  ambi¬ 
ent  magnetic  field  compensation  may  require  a  reference  magnetometer  so  close  to 
the  test  magnetometer  that  it  may  be  influenced,  too,  by  the  ship’s  magnetic  field. 
We  term  the  technique  applied  in  this  case  the  Local  Reference  Technique.  Fig.  2 
(red  lines)  demonstrates  its  application.  OEMs  03  and  04  are  both  influenced  by 
the  magnetic  field  of  Manning,  the  residual  between  the  two  OEMs  is  therefore  not 
identical  to  the  true  magnetic  field  perturbation,  specifically,  it  is  of  smaller  ampli¬ 
tude.  The  residual  does  thus  not  provide  for  a  true  replica  of  Manning’s  magnetic 
signature.  It  is  one  of  the  objectives  of  this  report  to  show  that  the  choice  of  the 
proper  reference  technique  depends  on  the  area  which  is  supposed  to  be  kept  under 
surveillance. 
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Data  Analysis  Method 


In  performing  the  data  analysis  we  proceeded  with  the  following  steps.  After  having 
checked  the  magnetic  field  and  housekeeping  data  firom  the  OEMs  for  consistency, 
the  recorded  time  stamps  were  validated.  The  clocks  of  the  OEMs  had  been  syn¬ 
chronized  with  a  GPS  receiver  before  deployment.  At  the  end  of  the  sea  trial,  after 
recovery  of  the  OEMs,  their  internal  clocks  were  compared  against  a  GPS  reference 
clock.  The  differences  between  nominal  and  actual  times  are  listed  in  Table  3. 


Table  3  OBM  oscillator  deviations  from  GPS  standard 


Year 

OBM 

Operation 

Clock  offset 

Normalized 

1996 

M03 

8  d 

-755  ms 

-94  ms/d 

M04 

20  d 

+125  ms 

+6  ms/d 

M06 

20  d 

-220  ms 

-11  ms/d 

M12 

20  d 

+1230  ms 

+62  ms/d 

M17 

16  d 

-75  ms 

-5  ms/d 

M18 

20  d 

-175  ms 

-9  ms/d 

M19 

20  d 

+65  ms 

+3  ms/d 

1997 

M04 

16  d 

-100  ms 

-6  ms/d 

M12 

16  d 

+590  ms 

+37  ms/d 

M17 

16  d 

-530  ms 

-33  ms/d 

M18 

16  d 

-500  ms 

-31  ms/d 

M19 

16  d 

-575  ms 

-36  ms/d 

All  housekeeping  data  were  exhaustively  scrutinized,  not  just  spotchecked.  It  was 
verified  that  the  time  stamps  were  always  correctly  updated,  the  temperatiure  re¬ 
mained  stable  and  at  a  meaningful  value,  the  supply  voltages  were  in  the  correct 
range  and  decreased  with  time  progressing  (though  very  little),  the  compass  bearing 
and  the  inclination  angles  remained  stable  once  the  modules  had  settled  on  the  sea 
bottom  (usually  within  a  few  hours  after  deployment),  and  the  magnetometer  x-  and 
y-ax.es  deviated  from  the  horizontal  plane  by  not  more  than  a  few  degrees.  Loca¬ 
tions  had  been  selected  for  deployment  where  the  sea  bed  was  practically  horizontal, 
and  this  was  verified  by  Manning’s  echo  sounder. 
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Once  the  raw  data  had  been  transferred  to  an  IBM  type  PC,  the  three  components 
of  the  magnetic  field  vector  were  transformed  into  a  suitable  coordinate  system. 
The  mathematical  formulation  of  the  rotation  procedure  is  given  in  Annex  A.  The 
1997  measurements  were  conducted  far  away  fi:om  the  mainland  where  no  implied 
preference  direction  of  the  magnetic  field  variations  exists.  We  adapted  a  coordinate 
system  aligned  with  the  geomagnetic  north,  east  and  vertical  axes.  For  the  near 
coast  measmements  of  1996,  we  found  that  the  amplitude  of  the  magnetic  field 
variations  was  related  to  the  orientation  of  the  coast  line,  with  mean  amplitudes 
different  between  the  components  parallel  and  perpendicular  to  the  coast.  Therefore 
it  appeared  appropriate  to  use  a  coordinate  system  in  which  the  x-axis  runs  parallel 
to  the  mean  coast  line,  the  j/-axis  perpendicular  to  it,  and  the  z-axis  points  vertically 
down. 

The  magnetic  field  data  were  low-pass  filtered  using  a  0.08  Hz  (12.5  s)  cutoff  fre¬ 
quency,  decimated  to  6  seconds  sampling  rate,  and  subsequently  high-pass  filtered 
with  a  cutoff  frequency  of  0.0015  Hz  (667  s).  This  bandwidth  was  chosen  such  that 
it  would  normally  encompass  the  magnetic  signature  of  a  passing  submarine  (c.f. 
Fig.  1).  Considering  that  the  clock  deviations  from  their  nominal  values  (Table  3) 
were  always  at  least  an  order  of  magnitude  smaller  than  the  low-pass  cutoff  period, 

no  need  arose  for  a  correction  of  the  OBM  time  stamps. 

)  ■■ 

The  magnetic  field  variations  observed  close  to  the  coast  are  largely  determined  by 
train  noise  from  a  busy  coastal  railway  line  (180  trains  per  day)  which  tends  to  follow 
a  pronotmced  diurnal  pattern.  Diurnal  magnetic  field  intensity  distributions  were  de¬ 
rived  from  16  days  of  continuous  land-based  measurements  near  Deiva  Marina,  made 
by  a  former  SACLANTCEN  summer  student,  A.  Magunia,  in  September/October 
1996,  during  our  sea  trial.  The  result  is  displayed  in  Fig.  5  in  a  self-normalized 
gray  scale  (A.  Magunia,  J.W.Goethe-Universitt,  Prankfurt/Medn,  private  commu¬ 
nication).  Based  on  this  distribution  of  the  magnetic  field  disturbance  level,  more 
specifically,  the  magnetic  field  variance  over  subsequent  15-min  sections,  each  day 
was  devided  into  four  intervals  which  are  listed  in  Table  4. 


Table  4  Diurnal  intervals  according  to  train  traffic 

Time  interval  (UTC) 

Characteristics 

04:30  -  09:30 

09:30  -  18:00 

18:00  -  21:30 

21:30  -  04:30 

moderate  train  traffic,  medium  noise 
heavy  train  traffic,  high  noise 
moderate  train  traffic,  medium  noise 

little  train  traffic,  low  noise 

i 


For  each  OBM,  each  day  of  operation,  and  each  of  the  intervals  listed  in  Table  4, 
the  recorded  magnetic  field  variation  was  centred  to  yield  zero  mean.  Subsequently, 
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the  amplitude  difference  between  two  neighbouring  OEMs  was  computed.  This 
method  is  a  simple  application  of  a  reference  technique.  It  was  improved  upon  in  a 
subsequent  step. 

Several  error  sources,  partly  device-dependent,  partly  imposed  by  nature,  render 
the  simple  vector  difference  an  unsatisfactory  method.  Firstly,  the  OEM  compass 
has  only  2°  accuracy,  i.e.  2°  uncertainty,  and  the  inclination  measurements  have  0.2° 
resolution  and  uncertainty.  Secondly,  the  three  axra  of  the  fluxgate  sensors  may  not 
be  exactly  orthogonal  and  aligned  with  the  compass  and  inclinometer  axes.  Thirdly, 
the  geography  (i.e.  the  distribution  of  water  and  land  masses)  and  the  geology  pro¬ 
duce  discontinuities  of  the  electrical  conductivity  distribution  and  impose  certain 
preferred  directions  on  those  electric  currents  which  are  induced  by  the  external 
magnetic  field  variations.  They  can  contribute  significantly  to  the  total  magnetic 
field.  Fourthly,  the  geometry  of  stray  currents  resulting  from  leakage  of  the  train 
return  current  is  largely  dependent  on  the  location  of  the  railway  power  substations, 
the  orientation  of  the  railway  line,  and  the  conductivity  structmre  of  the  seabed,  and 
may  thus  have  different  preference  directions  at  different  OEM  locations. 

In  order  to  alleviate  the  effect  of  these  errors,  a  time-domain  trivariate  linear  re¬ 
gression  algorithm  was  applied  to  the  vector  time  series  from  pairs  of  neighbouring 
OEMs  for  each  of  the  fom:  diurnal  intervals  and  each  day  of  magnetometer  operation. 
Such  a  method  had  also  been  applied  to  land-based  magnetic  field  measurements  in 
order  to  determine  the  impact  of  various  error  sources  on  the  apparent  spatial  cor¬ 
relation  between  two  adjacent  sites  [11].  The  regression  scheme  relates  the  magnetic 
field  variation,  observed  by  a  tri-axial  sensor  at  a  test  site,  to  the  magnetic  field 
variation,  observed  at  the  same  time  by  a  tri-axial  sensor  at  a  reference  site 
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The  residual  vector,  <JB^,  is  just  of  that  paxt  of  the  magnetic  field  that  can  not  be 
explained  through  linear  regression  between  test  and  reference  sensor.  If  no  instru¬ 
ment  error  or  magnetic  field  non-uniformity  would  affect  the  measurements  and  if 
the  three  magnetic  field  vector  components  were  fully  independent  from  each  other, 
{Aij)  were  the  identity  matrix.  Increasing  deviation  of  (Ay)  from  a  diagonal  matrix 
is  equivalent  to  an  increasing  cross-correlation  between  different  vector  components. 
The  coefficients  of  the  regression  matrix,  (Ay),  are  found  by  minimizing  the  sum 
of  the  squared  error  residuals  (the  elements  of  the  right-most  vector  on  the  right 
hand  side  of  eq.  (4)).  This  procedure  results  in  magnetic  field  vectors  which  we 
call  "optimally  aligned  and  scaled” .  The  mathematical  procedure  for  computing  the 
coefficients  is  described  in  Annex  B. 
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If  the  magnetic  field  variations  observed  simultaneously  by  two  magnetometers  con¬ 
sist  of  a  superposition  of  signals  highly  correlated  between  the  two  sites,  and  ran¬ 
dom  noise  of  equal  variance,  uncorrelated  between  the  two  sites,  the  error  terms  are 
smaller  if  the  magnetometer  with  the  more  intense  signal  is  taken  as  the  reference 
sensor.  This  is  shown  in  Annex  C.  In  the  vicinity  of  the  coast,  the  train  noise  is 
the  dominating  magnetic  signal,  the  intensity  of  which  decreases  with  increasing 
distance  from  coast.  Fig.  6  shows  a  stack  of  sample  time  series  firom  a  3-hour  inter¬ 
val,  recorded  in  the  fall  of  1996  by  the  OEMs  17,  18,  19  and  12  (c.f.  Map  1).  Only 
the  coast-perpendicular  component  is  displayed.  We  notice  that  the  magnetic  field 
variations  axer  very  similax  in  shape,  and  the  amplitudes  decrease  with  increasing 
distance  from  the  coast.  Consequently,  we  chose  for  each  pair  of  magnetometers  the 
one  closer  to  coast  as  the  reference  sensor.  The  1997  measurements  were  made  far 
away  from  the  coast,  with  a  group  of  magnetometers  close  to  Montecristo  island  and 
another  group  off-shore  several  times  the  extent  of  the  island.  We  use  OEM  19,  one 
of  the  latter  group,  as  a  reference.  The  pairs  of  OEMs  selected  for  analysis  are  listed 
in  Table  5  (see  also  Maps  1  and  2).  The  meaning  of  Lobs  (last  column  of  Table  5) 
is  explained  in  Annex  D. 


Table  5  Pairs  of  OEMs  seltcted  for  regression  analysis 


Year 

Test->  Reference 

Distance 

Connection  line 

f^obs 

1996 

M18->^M17 

1.55  km 

perpendicular  to  coast 

0.4825  nT 

M19^M18 

1.97  km 

perpendicular  to  coast 

0.5455  nT 

M12^.M19 

2.40  km 

perpendicular  to  coast 

0.6100  nT 

M06->M12 

5.97  km 

parallel  to  coast 

1.1455  nT 

M04->M06 

7.94  km 

parallel  to  coast 

1.4410  nT 

M03^M12 

18.73  km 

parallel  to  coast 

1997 

M17->M19 

1.11  km 

so^th  to  north 

M04-4M19 

12.30  km 

south-east  to  north-west 

For  each  pairs  of  neighbouring  magnetometers,  each  day  of  operation,  and  each 
diurnal  interval,  the  sums  of  the  squared  absolute  and  relative  error  residuals  were 
examined.  Time  segments  having  abnormally  large  residuals  were  not  used  in  the 
following  step,  namely  averaging  of  the  matrix  coefficients.  The  number  of  intervals 
discarded  was  very  small,  they  are  listed  in  Table  6.  The  matrix  coefficients  of 
the  remaining,  ’’clean”  segments,  were  average4  over  consecutive  days  so  that  one 
single  transfer  matrix  resulted  for  each  of  the  fpur  diurnal  intervals  and  each  pair 
of  OEMs.  These  averaged  matrices  were  used  in  constructing  univariate  transfer 
function  models  in  the  time  domain,  of  the  type  described  by  eq.  (4). 
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Table  6  Time  segments  excluded  from  regression  matrix  averaging 


Year 

Test^Reference 

Time  interval  (UTC) 

1996 

M18->M17 

Oct  08.  09:30-18:00 

M19-4M18 

Oct  06,  09:30-18:00 

M12^M19 

Oct  06,  09:30-18:00 

M12^M19 

Oct  10.  09:30-18:00 

M06^M12 

Oct  02,  21:30-Oct  03.  04:30 

M04^M06 

Oct  04.  18:00-21:30 

M04->M06 

Oct  09.  18:00-21:30 

M03->M12 

Oct  10,  18:00-21:30 

1997 

M17^M19 

July  23,  04:30-09:30 

M17^M19 

July  23,  09:30-18:00 

This  method  was  chosen  in  order  to  simulate  the  data  processing  in  a  ’’real”  subma¬ 
rine  detection  experiment.  In  such  an  experiment,  a  reliable  transfer  matrix  would 
be  determined  from  "clean”  time  segments,  i.e.  time  segments  not  contaminated  by 
passing  ships  or  certain  types  of  sensor  noise.  The  transfer  matrix  would  subse¬ 
quently  be  applied  in  real  time  to  magnetic  field  measurements  coming  in  from  two 
magnetometers,  in  order  to  identify  local  magnetic  field  perturbations  which  appear 
at  only  one  of  the  two  OBM  sites,  or  at  both  sites  but  with  different  amplitudes  and 
orientations. 

In  the  following  discussion  we  consider  one  of  the  two  OEMs  of  each  pair  the  ’’test 
magnetometer”,  measurements  from  which  are  confronted  to  those  from  a  ’’reference 
magnetometer”.  The  time  series  of  magnetic  field  variations  measured  by  the  test 
magnetometer,  their  difference  to  the  variations  measured  simultaneously  by  the 
reference  magnetometer,  and  the  residual  time  series  which  remain  after  trivariate 
representation  through  averaged  transfer  matrices,  were  binned  according  to  the 
scheme  described  below.  Different  bin  sizes,  0.2  nT  for  the  1996  data  and  0.05  nT 
for  the  1997  data,  were  chosen  in  order  to  account  for  the  largely  different  mean 
intensity  levels  of  the  magnetic  variations  recorded  in  the  two  different  areas. 

The  binning  was  done  individually  for  each  pair  of  OEMs,  each  vector  component, 
each  day  and  each  diurnal  interval  individually.  It  was  based  on  the  deviation  of 
the  data  points  from  the  zero  line,  irrespective  of  their  sign,  whether  positive  or 
negative.  This  means,  in  effect,  coimting  the  number  of  times  the  absolute  value 
increased  past  one  of  a  series  of  regularly  spaced  thresholds,  and  tally  the  results. 
Crossings  of  the  zero  level  were  not  counted.  Fig.  7  illustrates  the  binning  procedure 
for  a  fictitious  data  set  and  four  bin  thresholds  on  each  side  of  the  zero  line. 
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A  complete  set  of  plots  from  all  pairs  of  magnetometers  listed  in  Table  5  appears 
in  Fig.  8a-h.  Magnetic  field  variation,  difference  to  the  neighbour  OBM,  and  the 
residual  computed  for  each  vector  component,  are  plotted  in  separate  panels.  The 
days  are  incremented  along  the  abscissa,  with  each  day  broken  into  fom:  columns 
which  represent  the  four  diurnal  intervals,  and  the  bins  along  the  ordinate.  The 
decadic  logarithm  of  the  number  of  counts  at  each  bin  threshold,  normalized  to  a 
common  time  unit  of  one  hour,  is  displayed  using  the  color  code  found  at  the  bottom 
of  each  page. 


Table  7  Identified  disturbances  removed  from  the  1996  data 


Module 

Date 

Time  (UTC) 

Description 

M03 

Oct  10 

20:30-20:50 

vessel 

M04 

Sep  30 

08:55-09:10 

Manning 

Oct  03 

19:20-19:30 

vessel 

Oct  04 

13:30-13:40 

Manning 

Oct  07 

08:25-08:30 

vessel 

Oct  08 

15:10-15:15 

Manning 

Oct  09 

19:40-19:55 

vessel 

M06 

Sep  30 

09:30-09:40 

Manning 

Oct  03 

00:45-01:05 

vessel 

Oct  04 

13:00-13:10 

Manning 

Oct  04 

19:05-19:15 

vessel 

Oct  08 

14:35-14:45 

Manning 

Oct  11 

14:55-15:05 

vessel 

M12 

Sep  30 

10:10-10:20 

Manning 

Oct  04 

12:35-12:40 

Manning 

Oct  08 

14:10-14:20 

Manning 

Oct  10 

10:15-10:25 

ve^el 

M17 

Oct  04 

12:05-12:15 

Manning 

Oct  08 

13:40-13:50 

Manning 

M18 

Sep  30 

15:45-15:55 

Manning 

Oct  04 

12:15-12:20 

Manning 

Oct  08 

13:50-14:00 

Manning 

M19 

Sep  30 

15:35-15:45 

Manning 

Oct  04 

12:20-12:30 

Manning 

Oct  06 

14:30-14:40 

vessel 

Oct  08 

14:00-14:10 

Manning 

Oct  13 

20:40-10:45 

ves^l 
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The  residual  panel  is  split  into  two  sections,  the  upper  one  showing  the  results  from 
all  data,  the  lower  one  the  results  after  the  data  had  been  cleaned  for  those  artificial 
magnetic  field  perturbations  which  could  be  identified  with  certainty.  Removed  in 
this  process  were  signatures  of  passing  vessels  (including  site  visits  by  Manning), 
acoustic  interrogation,  individual  instrument  spikes  lasting  for  one  single  data  point 
only,  and  obvious  sensor  noise  of  large  amphtude,  such  as  abrupt  offset  changes. 
Tables  7  and  8  list  the  data  which  were  removed  fi:om  the  residual  time  series  in 
order  to  produce  the  cleaned  residuals. 


Table  8  Identified  disturbances  removed  from  the  1997  data 


Module 

Date 

Time  (UTC) 

Description 

M03 

July  23 

04:30-04:36 

07:24-07:42 

15:56-15:57 

18:55-19:05 

sensor  noise 

sensor  noise 
single  spike 
sensor  noise 

M04 

July  21 

14:36-14:40 

Manning 

M17 

July  20 

18:09-18:19 

vessel 

July  24 

09:40-09:50 

vessel 

M18 

July  24 

09:40-09:50 

vessel 

July  25 

04:35-04:45 

sensor  noise 

Mig 

July  21 

16:41-16:43 

17:07-17:08 

17:28-17:29 

single  spike 
single  spike 
single  spike 

July  23 

12:54-12:55 

17:14-17:16 

single  spike 
single  spike 
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6 

Distributions  of  Magnetic  Field 
Variations,  Differences  and  Residuals 


The  data  from  the  1997  cruise,  taken  in  an  area  far  away  from  anthropogenic  mag¬ 
netic  field  sources,  may  serve  as  a  manifestation  of  the  "ground  state”  of  magnetic 
field  fiuctuations,  against  which  the  data  taken  at  magnetically  highly  disturbed 
sites  along  the  Ligurian  coast  (1996  cruise)  have  to  be  confronted.  We  first  discuss 
two  data  sets  from  the  1997  cruise  (Fig.  8a,b)  and  then  a  series  of  data  sets  from 
the  1996  cruise  (Fig.  8c-h). 

Fig.  8a  shows  the  magnetic  field  variations  observed  by  OBM  17,  and  a  comparison 
with  the  reference  magnetometer  OBM  19.  In  the  X-  and  F-components  (north¬ 
ward  and  eastward)  we  notice  a  slight  dependence  of  the  magnetic  field  intensity 
on  the  time  of  the  day.  It  appears  to  be  highest  in  the  first  two  diurnal  intervals 
(04:30-18:00  UTC,  equivalent  to  06:30-20:00  local  time)  and  lowest  in  the  forth 
(21:30-04:30  UTC,  23:30-06:30  local  time).  Anthropogehic  electromagnetic  noise 
originating  from  Montecristo  island,  more  than  12  km  away  and  with  no  electro¬ 
magnetic  source  but  a  small  Diesel  generator,  can  not  be  responsible  for  the  diurnal 
pattern.  We  conclude  that  even  areas  of  the  Tyrrhenian  Sea  some  50  km  away 
from  the  mainland  are  to  a  small  extent  contaminated  by  anthropogenic  magnetic 
noise.  We  also  notice  that  the  vertical  magnetic  field  component  is  generally  much 
weaker  than  the  horizontal  one.  This  is  typical  for  a  horizontally  layered  electrical 
conductivity  distribution  with  a  good  conductor  on  the  top,  such  as  sea  water  and 
water-saturated  sediments.  The  magnetic  field  difierence  between  OBMs  17  and  19 
is  much  smaller  than  their  magnetic  field  variations  which  demonstrates  that  the 
variations  are  to  a  high  degree  correlated  between  both  sites.  The  even  smaller 
residual  amplitudes  confirm  this  but  uncover  also  several  exceptions. 

Two  segments  bearing  large  residuals,  the  second  diurnal  interval  of  July  20  and 
the  third  of  July  21,  were  found  to  be  affected  by  the  signature  of  a  vessel  passing 
OBM  17,  and  a  group  of  single  spikes  was  detected  in  the  data  of  OBM  19  (c.f. 
Table  8).  Otherwise,  the  residuals  remain,  with  one  exception  which  occured  during 
the  second  interval  of  July  20,  below  the  dual-instrument  maximum  noise  level 
(according  to  our  estimate  about  0.25  nT  over  the  frequency  band  considered  here). 
On  July  23,  the  data  from  OBM  17  appeared  highly  disturbed,  with  a  sequence 
of  magnetic  field  perturbations  bearing  all  the  characteristics  of  instrument  noise. 
Altough  the  magnetometer  tended  to  recover  somewhat,  it  never  returned  to  the 

SaCLANTCEN  SR-304 


-18- 


SaCLANTCEN  SR-304 


previous,  low  noise  level,  a  fact  which  is  particularly  well  seen  in  the  data  from  the 
Z-component  (vertical).  The  effect  is  well  visible  in  the  difference  and  residuals,  in 
fact,  from  July  23  on,  the  Z-component  difference  and  residual  remain  almost  as 
large  as  the  Z  variation  of  OBM  17.  It  means  that,  beginning  with  July  23,  the 
correlation  between  the  data  from  the  two  magnetometers  became  very  small.  A 
possible  but  not  confirmed  explanation  could  be  the  following.  The  modules  were 
exposed  to  direct  sunlight  and  heated  up  while  being  prepared  for  deployment  from 
the  workboat,  and  a  valve  had  to  be  opened  to  release  the  air  pressure  inside  the 
module  which  started  to  build  up  with  the  rise  of  the  temperature  (36°  C  peak 
measured  immediately  before  deployment).  The  warm,  humid  air  above  the  water 
had  thus  access  to  the  interior  of  the  modules.  Once  deployed  and  adjusted  to 
the  significantly  lower  temperature  prevailing  on  the  sea  bottom  (16°  C  measured), 
the  humid  air  inside  may  have  condensed  to  drops  which  affected  negatively  the 
electronic  circuitry,  despite  precautions  taken  against  humidity. 

Fig.  8b  shows  further  results  from  the  same  cruise.  Data  from  OBM  04  (located 
close  to  Montecristo  island)  are  compared  to  those  from  OBM  19,  i.e.  over  a  dis¬ 
tance  of  more  than  12  km.  We  notice  that  the  residuals  between  both  data  sets,  after 
removing  the  signature  of  Manning  and  the  sequence  of  instrument  spikes  from  the 
second  interval  of  July  21,  rarely  exceed  0.5  dT.  A  comparison  with  the  first  four 
undisturbed  days  of  the  OBM  pair  17/19  (Fig.  8a)  shows  that  the  residuals  OBM 
04  vs.  OBM  19  eire  only  twice  as  large  as  the  residuals  OBM  17  vs.  OBM  19.  Only 
around  July  24,  intense  magnetic  activity  of  ionospheric  origin  led  to  large  natural 
magnetic  field  variations  which  left  their  traces  on  the  residuals,  making  them  larger 
than  they  were  during  magnetically  quiet  times.  We  suggest  that  the  area  investi¬ 
gated  in  1997  is  suitable  for  application  of  the  Remote  Reference  Technique,  with  a 
sensor  spacing  of  the  order  of  2  to  3  kilometres. 

We  now  turn  to  the  data  from  the  1996  Ligurian  coast  cruise.  Fig.  8c-h.  They  reveal 
various  characteristics  which  are  common  to  all  of  the  OBMs. 

•  A  pronounced  and  regularly  repeated  diurnal  pattern  is  observed  in  all  three 
magnetic  field  components.  Amplitudes  of  magnetic  field  variations  recorded 
throughout  the  1996  cruise  are  much  larger  than  those  from  the  1997  cruise. 

•  The  Y -component  (oriented  coast-perpendicular)  generally  exhibits  the  small¬ 
est  magnetic  field  variations,  and  the  Z-component  (oriented  vertical)  the 
largest. 

•  When  inspecting  magnetic  field  vector  differences  and  residuals  between  two 
neigbouring  OBMs,  the  Z-component  remains  to  be  the  largest  one. 

The  following  properties  apply  only  to  the  OBMs  04,  06,  12,  17,  18  and  19.  OBM 
03,  much  farther  to  the  northwest,  has  to  be  exempted. 
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•  The  vector  differences  between  neighbouring  magnetometers  (middle  row  pan¬ 
els  of  Fig.  8c-g)  are  smaller  than  the  variations  (top  row  panels)  and  facilitate 
detecting  signatures  of  local  magnetic  sources  which  affect  only  one  of  the  two 
OBMs  in  a  pair. 

•  The  residuals  seem  to  improve  upon  the  differences  in  the  sense  that  they 
are  significantly  smaller  than  the  differences,  in  particular  in  the  X-  and  Y- 
components,  less  so  in  the  .^-component.  This  facilitates  even  more  detecting 
local  magnetic  field  perturbations. 


OBM  03  differs  firom  the  rest  of  the  array  in  a  particular  way: 


•  The  magnetic  field  differences  to  the  magnetometer  closest  by,  OBM  12,  ex¬ 
ceed  the  magnetic  field  variations,  and  the  residujals  are  of  the  same  order 
of  magnitude  as  the  variations.  Magnetic  field  variations  at  the  sites  of  the 
OBMs  03  and  12  appear  to  be  poorly  correlated. 


The  following  propositions  are  made  in  order  to  explain,  at  least  partially,  the  find¬ 
ings  described  above. 

Natural  magnetic  field  fluctuations  are  solar-wind  driven  and  are  unevenly  and  irreg¬ 
ularly  distributed  over  the  day.  At  midlatitude  stations,  narrow-banded,  continuous 
types  of  pulsations  prevail  on  the  sunlit  side  of  the  earth,  and  irregular,  bursty  types 
on  the  dark  side.  Our  findings  from  our  coastal  OBM  array  suggest,  however,  a  dom¬ 
inating  influence  of  anthropogenic  sources  with  a  rather  regular  diurnal  variation. 
The  Italian  state  railway  system,  which  uses  DC-powered  electric  engines,  is  the 
most  likely  candidate  for  their  generation  [12],  [13],  [14]. 

The  train  engines  are  powered  from  substations  located  along  the  railway  line  at 
a  spacing  of  10  to  30  kilometres.  Part  of  the  retiurn  current  from  the  trains  flows 
through  the  rails  back  to  the  substations,  but,  due  to  intentional  grounding  of  the 
rails  at  regular  spacing,  another  part  of  the  retm^n  current  leaks  and  finds  other 
paths.  Each  time  a  train  accelerates  or  decelerates,  the  electric  power  consumption 
changes  drastically,  and  the  return  cmrent  to  the  power  substation  changes  with  it. 
The  high  intensity  of  magnetic  field  variations  observed  several  kilometres  from  the 
coast  (and  from  the  railway  line)  indicates  that  leakage  currents  reaching  far  out 
into  the  sea  must  be  the  source  of  the  major  part  of  the  magnetic  field  variations.  A 
simple  calculation  reveals  that  the  electric  ciurent  loop  formed  by  the  supply  current 
through  the  overhead  cont2«;t  wire,  the  train  engine,  and  the  rail  return  current,  can 
not  account  for  measurable  magnetic  field  variations.  Using  t)rpical  numbers,  namely 
a  supply  cmrrent  of  1000  A  and  a  vertical  distance  between  contact  wire  and  rail  of 
5  m,  we  find  that  the  magnetic  signal  strength  2  km  away  (i.e.  at  the  location  of  the 
nearest  module,  OBM  17)  should  be  0.25  nT,  equal  to  the  dual-sensor  self  noise. 
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8  km  away  from  the  railway  (i.e.  at  the  location  of  OBM  12),  the  signal  would  be 
only  15  pT  and  thus  undetectable.  We  conclude  that  the  leakage  currents  must  flow 
much  closer  to  the  OBMs  than  the  rails  run  (i.e.  the  ciurrents  extend  far  out  into  the 
sea)  in  order  to  generate  the  large  magnetic  fleld  variations  observed  by  the  OBM 
array. 

Let  us  now  consider  a  northwestbound  train  starting  from  La  Spezia.  Before  reaching 
Levanto  the  train  is  powered  by  the  substations  of  La  Spezia  and  Levanto.  La 
Spezia  is  located  on  the  remote  side  of  the  mountain  ridge,  and  the  electrically 
poorly  conducting  mountains  largely  prevent  leakage  currents  between  La  Spezia 
and  Levanto  to  reach  the  sea  where  the  OBM  array  resides.  Thus,  only  leakage 
currents  between  Levanto  and  the  train  play  a  role  for  OBM  04.  When  running 
between  Sestri  Levante  and  Levanto  the  train  engine  is  powered  by  the  substations 
at  the  latter  two  cities,  and  the  return  leakage  currents  flow  back  to  these  two 
substations.  OBMs  17,  18,  19,  12  and  06  are  thus  simultaneously  exposed  to  the 
cmrrents  supplied  by  and  flowing  back  to  Levanto  and  Sestri  Levante,  but  OBM  03 
is  not.  Once  the  train  passes  Sestri  Levante,  the  train  engine  draws  power  from 
the  substations  of  Sestri  Levante  and  San  Margherita  Ligure,  and  now  OBM  03  is 
exposed  to  train  leakage  currents  while  the  rest  of  the  array  is  no  longer  affected 
by  it.  Though  the  large  magnetic  field  differences  and  residuals  between  OBM  03 
and  the  rest  of  the  array  may  be  explained  this  way,  we  still  expect  a  noticable 
difference  between  OBMs  04  and  06.  This  is  not  the  case,  the  pairs  OBM  04/06 
and  OBM  06/12  have  rather  similar  difference  and  residual  levels,  both  of  which  are 
much  smaller  than  the  variations. 

Fig.  9  gives  an  example  of  the  magnetic  disturbance  caused  by  the  passage  of  a  single 
train.  During  the  time  interval  shown  here,  only  one  single  train  ran  between  Levanto 
and  Sestri  Levante,  according  to  the  records  of  the  regional  railway  administration 
in  Genoa  (we  did  not  rely  on  the  official  timetable  which  may  be  inaccurate).  The 
highlighted  sections  show  when  the  train  was  on  the  rails  between  Levanto  and  Sestri 
Levante  (yellow)  and  when  it  traversed  the  railway  station  of  Deiva  Marina  (green). 
The  train  did  not  stop  there,  but  the  large  magnetic  field  offsets  indicate  that  the 
train  did  significantly  decelerate  before  and  accelerate  after  the  passage. 

The  coastal  effect,  which  explains  the  prevention  of  large-scale  electric  currents, 
such  as  currents  induced  by  natural  magnetic  fields,  from  penetrating  from  the  well 
conducting  sea  into  the  poorly  conducting  coastal  mountain  range,  would  result  in 
large  vertical  magnetic  field  variations  close  to  the  coast  [15].  When  comparing  the 
upper  right  panels  of  Fig.  8c-g  one  notices  some  attenuation  of  the  Z  variations  with 
increasing  distance  from  the  coast  (Fig.  8c-e)  and  practically  equal  amplitudes  at 
sites  having  about  the  same  distance  from  coast  (Fig.  8e-g).  This  observation  is 
consistent  with  the  coastal  effect. 
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The  horizontal  magnetic  field  variations  show  large  amplitudes  in  the  X-component 
(coast-parallel),  which  decrease  with  increasing  distance  firom  coast  and  remain  ap¬ 
proximately  constant  at  constant  distance  from  coast  (c.f  Fig.  8c-g,  upper  left  pan¬ 
els),  and  less  intense  variations  in  the  F-component  (coast-perpendicular).  Electric 
leakage  currents  from  the  railway  line,  if  fiowing  close  to  the  rails,  would  produce 
a  magnetic  field  mainly  perpendicular  to  the  coast  (i.e.  in  the  Y-  rather  than  the 
X-component).  Our  observations  thus  imply  that  the  railway  leakage  currents  reach 
far  out  into  the  sea  and  confirm  the  statement  we  made  above. 
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7 

Cumulative  Binning  of  the  Distributions 


The  binned  data  which  had  been  used  to  create  Fig.  8a-h,  were  then  processed  in 
a  different  way.  Instead  of  creating  distributions  for  each  day  and  diurnal  interval 
individually,  all  data  from  a  given  pair  of  magnetometers  and  a  given  diurnal  in¬ 
terval  were  accumulated  and  then  normalized  to  yield  threshold  crossings  per  hour 
(’’events/hour”).  Only  the  OBM  pair  17/19  (1997  sea  trial)  was  divided  into  two 
sections,  one  prior  to  July  23,  and  the  other  starting  from  July  23.  This  was  deemed 
necessary  because  of  the  significantly  different  behaviom  of  OBM  17  from  July  23 
on  (see  discussion  above).  The  OBM  pair  03  12  (1996  sea  trial)  is  not  included  in 
Fig.  11  because  Fig.  8h  demonstrates  that  the  residuals  were  as  large  as  the  varia¬ 
tions  themselves,  i.e.  there  is  very  little  correlation  between  the  two  sites.  No  further 
comparison  is  deemed  useful. 

To  facilitate  inter-comparison  of  different  OBM  pairs  from  1996,  we  normalized 
the  binning  results  to  a  distance  of  1  km,  assuming  a  linearly  increasing  difference 
and  residual  magnetic  field  and  constant  sensor  noise.  In  other  words,  we  display 
averages  of  the  spatial  gradients  of  the  magnetic  field  differences  and  residuals,  sup¬ 
plemented  by  system  noise.  In  essence,  the  distributions  were  modified  according  to 
Annex  D,  eq.  (D3).  This  procedure  incorrectly  modifies  the  magnetic  field  variations 
which  are  per  definition  single  station  measures  and  for  which  the  distance  to  the 
neighbouring  OBM  bears  no  meaning.  In  particular,  a  division  of  the  amplitude 
by  the  distance  is  meaningless.  However,  we  accept  this  incorrectness  in  Fig.  11 
because  the  differences  and  residuals  rather  than  the  variations  are  the  important 
parameters  for  our  discussion. 

The  results  of  the  ciunulative  binning  are  shown  in  Fig.  10  and  11.  The  residuals 
used  in  generating  Fig.  10  and  11  are  the  ’’cleaned”  ones,  obtained  after  removing  the 
pertmbations  listed  in  Tables  7  and  8  (see  also  discussion  of  Fig.  8a-h).  Fig.  lOa-c 
refers  to  the  1997  (Montecristo)  sea  trial.  Fig.  lla-e  to  the  1996  (Ligurian  coast)  sea 
trial.  Each  page  shows  data  from  a  peirticular  pair  of  OBMs.  Each  vector  component 
and  each  diurnal  interval  occupies  one  panel  on  the  page.  The  number  of  events  per 
homr  (more  precisely,  their  decadic  logarithm)  is  displayed  using  different  colors  for 
the  variations,  the  differences,  and  the  residuals. 

Fig.  10  and  11  may  be  used  in  the  following  way:  For  a  selected  pair  of  magne¬ 
tometers  (which  means  basically  a  specific  geographic  location),  one  refers  to  the 
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corresponding  page.  For  a  selected  diurnal  interval  (out  of  four  possible)  and  a  se¬ 
lected  magnetic  field  component  (X,  Y  or  Z),  one  looks  up  the  corresponding  panel. 
For  a  selected  sensitivity  threshold  (in  units  of  nT  and  nT/km,  respectively),  one 
finds  in  the  panel  the  log-scaled  firequency  (in  events  per  hour)  of  the  magnetic  field 
differences  and  residuals  exceeding  the  specified  threshold. 

Qualitatively,  nothing  new  is  learned  fi"om  Fig.  10  and  11  when  confronted  with 
Fig.  8.  Quantitatively,  the  following  type  of  statements  can  be  made.  The  upper 
left  panel  of  Fig.  10a  may  serve  as  a  first  example.  Let  us  assume  that  magnetic 
field  measmrements  were  made  every  day  but  only  between  04:30  and  09:30  UTC, 
with  two  properly  operating  OEMs,  17  and  19.  We  find  tjie  decadic  logaxithm  of  the 
number  of  events/hom:  of  the  X  (north)  component  of  the  magnetic  field  variation  of 
OEM  17  to  be  about  -0.7  at  1.3  nT  bin  size,  -1-0.15  at  0.8  nT  bin  size,  and  as  much 
as  -1-1.05  at  0.4  nT  bin  size.  This  component  could  thus  be  expected  to  pass  once 
every  5  hours  a  level  of  ±1.3  nT,  once  every  40  minutes  a  level  of  ±0.8  nT,  and  about 
every  5-6  minutes  a  level  of  ±0.4  nT.  The  logarithm  of  the  number  of  events/hour 
of  the  X-component  of  the  magnetic  field  difference  between  OEMs  17  and  19  is 
-1.2,  i.e.  it  would  exceed  ±0.4  nT  only  once  dmring  16  hours  of  observations.  The 
X-component  magnetic  field  residual  would  practically  never  exceed  ±0.25  nT,  i.e. 
it  is  basically  the  magnetometer  self  noise  which  sets  the  limit  for  the  residuals. 

Fig.  10a  can  be  interpreted  in  the  following  way.  If  the  operator  of  an  OEM  array 
residing  in  the  vicinity  of  Montecristo  island  would  set  a  detection  threshold  of 
±0.4  nT,  the  X-component  of  a  single  tri-axial  magnetometer,  even  if  properly 
functioning,  would  issue  some  55  false  alarms  every  day  between  04:30  and  09:30 
UTC.  If  the  X-component  difference  with  a  reference  magnetometer,  1.1  km  away, 
were  considered,  the  operator  would  receive  a  false  alarm  every  three  days  at  most, 
and  if  the  operator  would  consider  the  residual  and  make  smre  that  he  has  become 
aware  of  all  passings  of  surface  vessels,  he  would  never  be  fooled  by  a  false  alarm. 
Such  a  statement  depends,  of  course,  on  the  area  and  the  electromagnetic  conditions 
prevailing  therein. 

To  demonstrate  this  we  take  Fig.  11a  for  a  second  example.  As  before,  we  consider 
only  mecisurements  made  during  the  04:30-09:30  UTC  interval  and  examine  only  the 
X-component.  All  three  components  of  the  magnetic  fi.eld  variation  of  OEM  18  reach 
several  nanoTesla  many  times  per  hour  and  are  thus  iiseless  for  MTE  purposes.  The 
difference  to  the  neighbouring  sensor,  OEM  17,  exceeds  ±0.74  nT/km  once  dtiring  4 
homrs,  and  ±0.48  nT/km  once  every  25  minutes.  The  residual  between  OEMs  18  and 
17  woiild  exceed  ±0.35  nT/km  once  every  12.5  hours  but  never  exceed  ±0.48  nT/km. 
If  the  operator  would  set  a  detection  threshold  oif  ±0.48  nT /km,  not  too  different 
from  the  one  used  above,  he  would  receive  a  false  alarni  more  than  twice  every  hour 
if  he  pays  attention  to  the  magnetic  field  difference,  but  he  would  practically  never 
receive  a  false  alarm  if  he  used  the  residual  instead.  The  succeeding  diurnal  interval, 
09:30-18:00  UTC,  is  much  more  disturbed.  If  we’cOMider  the  same  OEM  pair  and 
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the  same  vector  component  we  notice  that  the  difference  exceeds  a  threshold  of 
±0.74  nT/km  every  100  minutes  and  a  threshold  of  ±0.48  nT/km  every  12  minutes. 
The  residual  exceeds  a  threshold  of  ±0.35  nT/km  once  every  4  hours  and  a  threshold 
of  ±0.48  nT/km  every  19  hours  (left  panel,  second  row  of  Fig.  11a). 

Further  characteristics  of  the  data  can  be  deduced  from  Fig.  11.  The  magnetic 
field  variations  axe  generally  larger  in  the  coast-parallel  component  than  in  the 
coast-perpendicular  one,  eind  much  larger  in  the  vertical  component  than  in  the 
coast-parallel  one.  The  differences  and  residuals,  however,  axe  generally  smaller  in 
the  coast-parallel  component  than  in  the  coast-perpendicular  one,  and  still  quite 
laxge  in  the  vertical  component.  A  magnetometer  axray  operator  would  thus  pay 
little  attention  to  an  alert  issued  by  the  vertical  component,  and  most  attention  to 
one  issued  by  the  coast-parallel  component. 

This  should,  however,  not  be  interpreted  as  an  endorsement  for  the  use  of  two-axis 
magnetometers.  Firstly,  the  magnetometer  resides  in  an  eirbitraxy  attitude  on  the 
sea  bottom,  and  all  three  components  are  needed  in  order  to  realign  the  vector 
measurements  with  a  common  coordinate  system,  such  eis  coast-parallel,  coast- 
perpendiculax,  vertical.  Secondly,  tests  with  the  magnetic  field  data  used  in  this 
report  revealed  that  the  residuals  of  the  horizontal  magnetic  field  components  be¬ 
come  smaller  if  the  vertical  component  is  included  in  the  scheme  (trivariate  linear 
regression)  as  opposed  to  a  scheme  which  relies  exclusively  on  the  horizontal  com¬ 
ponents  (bivariate  linear  regression). 
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8 

Conclusion 


We  have  demonstrated  that  a  local,  time-dependent  perturbation  of  the  ambient 
magnetic  field,  as  it  is  produced  by  a  moving  submarine,  h^  to  compete  with  mag¬ 
netic  field  fluctuations  firom  other,  natural  and  anthropogenic,  soiurces.  The  latter 
are  often  much  more  intense  than  the  signal  from  a  passing  submarine.  The  cru¬ 
cial  parameter  which  determines  the  outcome  of  the  competition  is  the  spatial  scale 
length  of  the  magnetic  field  pertmrbations  in  the  bandwidth  under  consideration. 

The  bandwidth  under  consideration  is  limited  by  the  minimum  and  maximum  time 
interval  a  submarine  would  t3rpically  need  to  traverse  the  sphere  of  influence  of  the 
magnetic  field  sensor.  This  depends  on  such  parameters  as  the  system  noise  of  the 
sensor,  the  size  of  the  submarine,  and  its  magnetic  treatment.  For  the  Ocean  Bottom 
Magnetometers  (OEMs)  used  at  SACLANTCEN,  and  a  medium-sized,  imtreated 
submarine,  the  radius  of  the  sphere  of  influence  would  be  of  the  order  of  600  metres 
in  the  bandwidth  used  here,  0.0015-0.08  Hz. 

If  the  spatial  scale  length  of  the  non-submarine  magnetic  field  is  larger  than  that  of 
the  submarine  magnetic  field,  a  reference  technique  can  be  applied  which  allows  to 
suppress  the  fluctuations  of  the  ambient  magnetic  field  by  comparison  with  spaced 
reference  sensors.  The  spatial  scale  length  depends  largely  on  the  geographic  area 
where  the  observations  are  made.  In  this  report  we  have  discussed  data  from  two 
different  areas,  one  magnetically  very  uniform  (the  Montecristo  sea  area),  the  other 
magnetically  highly  non-uniform  (the  Ligurian  coastal  water). 

The  Montecristo  area  can  be  dealt  with  in  a  rather  simple  way;  the  spatial  correlation 
between  two  OEMs  is  only  limited  by  the  noise  characteristics  of  the  sensor  pair, 
even  over  a  distance  of  several  kilometres,  by  far  exceeding  the  submarine  sphere 
of  influence.  Fig.  8a  and  8b  demonstrate  this  by  showing  that  the  residual  (i.e. 
the  part  of  the  signal  that  is  not  linearly  correlated  between  two  seizors),  is  of 
the  order  of  0.25  nT,  virtually  independent  of  the  spacing  between  the  OEMs.  In 
this  area,  OEMs  could  be  deployed  with  a  ca.  1  km  spacing  which  would  gucurantee 
that  not  more  than  one  OEM  at  a  time  is  subject  to  the  submarine  field  while  all 
sensors  observe  practically  the  same  ambient  field  which  can  be  removed  by  linear 
regression  or  correlation  techniques.  The  Montecristo  area  is  therefore  a  suitable 
candidate  area  for  application  of  the  Remote  Reference  Technique. 

SaCLANTCEN  SR-304 


-26- 


SaCLANTCEN  SR-304 


The  Ligurian  coastal  water  is  subject  to  very  intense  and  localized  electromagnetic 
noise,  predominantly  from  a  busy  coastal  railway  line.  In  this  area,  the  spatial 
uniformity  of  the  ambient  magnetic  field  fluctuations  is  much  poorer  than  in  the 
Montecristo  area.  To  quantify  this,  we  formed  pairs  of  neighbouring  OEMs  and 
computed  the  magnetic  field  vector  differences  and  residuals.  The  residual  is  that 
part  of  the  magnetic  field  fluctuation  which  is  not  linearly  correlated  between  the  two 
OEMs  and  which  is  therefore  not  cancelled  out  when  applying  a  reference  technique. 

To  better  assess  the  diurnal  variation  ob  the  anthropogenic  magnetic  noise,  we  di¬ 
vided  the  day  into  four  diurnal  intervals  with  significantly  different  railway  activity 
(Table  4).  We  removed  all  clearly  identifiable  signatures  from  passing  surface  ves¬ 
sels  and  all  obvious  sensor  spikes  and  offset  jumps  from  the  data  and  computed 
magnetic  field  residuals,  using  one  specific  fixed  regression  matrix  for  each  pair  of 
magnetometers  and  each  diurnal  interval.  The  regression  matrices  had  been  deter¬ 
mined  from  a  correlation  analysis,  using  the  entire  data  set  available,  but  cleaned 
for  identified  ship  signatures  and  anomalous  sensor  noise  such  as  mentioned  above. 
We  then  determined  how  often,  for  each  pair  of  OEMs,  each  vector  component,  and 
each  diurnal  interval,  the  magnetic  field  variation,  difference,  and  residual  would 
exceed  a  given  threshold.  The  results  axe  graphically  displayed  in  Fig.  8. 

For  each  pair  of  neighbouring  OEMs  in  operation  during  the  Ligurian  coast  sea 
trial  and  each  magnetic  field  vector  component,  the  residual  was  normalized  to 
a  hypothetic  spacing  of  1  km,  taking  into  account  the  non-normalizable  system 
noise  contribution.  In  other  words,  we  computed  the  mean  spatial  gradient  of  the 
non-correlated  paxt  of  the  magnetic  field  fluctuations.  We  then  determined  the 
number  of  times,  averaged  over  the  entire  sea  trial,  which  this  mean  spatial  gradient 
would  exceed  a  given  threshold  (expressed  in  nT/km).  A  complete  summary  of  the 
results  is  found  in  Fig.  11.  In  order  to  clearly  distinguish  a  passing  submarine  from 
ambient  magnetic  noise,  the  occurence  frequency  of  suspicious  residuals  must  be 
made  low,  either  by  setting  the  threshold  high,  or  by  reducing  the  spacing  between 
magnetometers.  An  example  which  shows  how  to  use  the  results  from  the  previous 
sections  is  given  in  the  following. 

If,  for  instance,  a  magnetometer  array  operator  would  set  a  threshold  for  the  residual 
of  0.4  nT/km,  our  results  would  provide  him  with  the  probability  of  how  often 
he  would  find  this  threshold  exceeded  at  each  location,  without  the  influence  of  a 
submarine  or  surface  vessel.  The  result,  averaged  over  all  days  and  diurnal  intervals 
of  the  entire  sea  trial,  is  displayed  in  Fig.  12.  The  figure  shows,  for  each  pair  of 
OEMs  (18  17,  19  18,  etc.),  the  number  of  hours  that  would  typically  pass  before 
the  residual  exceeds  the  selected  threshold.  In  other  words,  if  the  operator  considers 
each  time  the  magnetic  field  residual  passes  the  threshold  a  submarine  detection 
event.  Fig.  12  provides  him  with  an  estimate  of  the  false  alarm  rate. 
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Fig.  12  reveals  some  characteristics  part  of  which  we  mentioned  already  in  a  previous 
section. 

(1)  The  vertical  component  of  the  magnetic  field  residual  exceeds  the  threshold  of 
0.4  nT/km  several  times  per  hour  and  can  thus  be  considered  irrelevant  for  subma¬ 
rine  detection  purposes. 

(2)  The  coast-perpendicular  component  of  the  magnetic  field  residual  becomes  a 
more  reliable  submarine  detection  parameter  the  farther  away  firom  the  coast  the 
OBM  pair  is  located. 

(3)  The  coast-perpendicular  component  of  the  magnetic  field  residual,  computed 
from  OEMs  having  the  same  distance  from  the  coast  (06^  12  and  04->06),  has  a 
very  low  false  alarm  rate.  Hundreds  to  thousands  of  horns  pass  before  the  threshold 
is  exceeded. 

(4)  The  coast-parallel  component  of  the  magnetic  field  residual  is  comparable  to  the 
coast-perpendicular  component  when  considering  pairs  with  OEMs  having  different 
distances  from  coast  (18 ->17,  19 ->18  and  12-)^  19),  but  not  when  considering  pairs 
having  the  same  distance  from  coast  (06— )^12  and  04-)‘06). 

An  operator  tasked  with  magnetic  surveillance  of  the  Ligurian  coastal  zone  would 
therefore  be  well  advised  to  heed  the  following  recommendations. 


•  The  vertical  magnetic  field  component  has  little  relevance  for  submarine  de¬ 
tection 

•  The  OEMs  should  be  deployed  in  such  a  way  that  pairs  can  be  formed  from 
OEMs  having  the  same  distance  from  coast. 

•  In  such  a  case,  the  coast-perpendicular  component  of  the  magnetic  field  resid¬ 
ual  yields  the  best  detection  parameter. 


We  emphasize  again  that  these  results  apply  only  to  the  particular  area  for  which  the 
study  was  made,  and  to  the  particular  magnetometers  used  for  the  measurements. 
Note  that  the  sensors  contribute  already  0.25  nT  to  the  hypothetical  threshold  of 
0.4  nT/km.  In  a  different  area,  the  results  may  be  different,  as  far  as  the  roles  of 
the  coast-perpendicular  and  coast-parallel  components  are  concerned.  Close  to  the 
coast,  the  vertical  component  will  always  be  a  poor  parameter,  for  two  reasons. 

(1)  Most  of  the  anthropogenic  electric  currents  flow  between  different  supply  points 
close  to  the  coast,  at  sea  surface  level,  i.e.  they  flow  horizontally  but  are  impeded  by  a 
lateral  electrical  discontinuity  (namely  the  coast).  It  is  mostly  the  vertical  magnetic 
field  which  reflects  strongly  the  horizontal  transition  from  a  good  conductor  to  poor 
conductor. 

(2)  The  coastal  effect  on  electromagnetic  induction  converts  temporal  variations 
of  the  horizontal  magnetic  field  into  vertical  variations.  As  far  as  magnetic  field 
correlation  is  concerned,  the  coastal  effect  also  favours  two  sites  having  the  same 
distance  from  coast  over  those  having  different  distances  from  coast. 
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Thus,  a  further  recommendation  to  the  OBM  array  operator  would  be: 


•  For  each  area  where  magnetic  surveillance  is  intended  to  be  established,  a  few 
days  of  magnetic  field  measurements  should  be  conducted,  prior  to  activating 
the  surveillance  system,  in  order  to  determine  which  vector  component  is  best 
suited  and  what  the  false  alarm  probability  is.  If  the  false  alarm  probability 
turns  out  to  be  unacceptably  high,  the  spacing  of  the  OEMs  should  be  de¬ 
creased  in  order  to  reduce  the  adverse  effect  of  a  large  spatial  gradient  of  the 
residual  magnetic  field. 


Decreasing  the  spacing  has  some  implication  on  the  reference  method  chosen.  Small 
sensor  spacing  (of  the  order  of  1  km  or  less)  requires  application  of  a  Local  Reference 
Technique^  large  spacing  permits  to  use  a  Remote  Reference  Technique. 

The  difference  between  the  two  techniques  is  demonstrated  in  Fig.  2  of  this  report 
and  is  here  briefly  recapitulated.  In  a  Remote  Reference  Technique,  only  one  of  the 
two  OEMs  forming  the  pair  is  influenced  by  the  submarine.  Therefore,  the  residual 
magnetic  field  contains  its  full  magnetic  signature.  In  a  Local  Reference  Technique, 
both  OEMs  of  the  pair  may  record  a  magnetic  signal  firom  the  submarine  (but 
usually  differently),  with  the  effect  that  the  signal  is  partially  correlated  between 
the  two  OEMs,  and  the  residuals  reveal  only  a  fraction  of  the  submarine’s  magnetic 
signature. 
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Map  1  The  southern  Ligurian  coastal  one  with  magnetometer  sites  marked  by  black  dots.  The  red  dot  in  the  inlay 
map  of  Italy  indicates  the  geographic  location  of  the  measurement  area.  Bathymetry  in  metres. 
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Map  2  Montecristo  island 
(Tyrrhenian  Sea)  and  its  easterly 
water  space,  with  magnetometer 
sites  marked  by  black  dots.  The  red 
dot  in  the  inlay  map  of  Italy 
indicates  the  geographic  location  of 
the  measurement  area.  Bathymetry 
in  metres. 


Map  3  The  Formiche  di  Grosseto 
area  (Tyrrhenian  Sea)  with  part  of 
Tuscany  visible  in  the  upper  right 
corner.  Magnetometer  sites  are 
marked  by  black  dots.  The  red  dot  in 
the  inlay  map  of  Italy  indicates  the 
geographic  location  of  the 
measurement  area.  Bathymetry  in 
metres. 
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Magnetic  field  distortion  functions 
(Anderson  functions) 


elapsed  time  (seconds)  relative  to  point  of  closest  approach 


(parameters:  speed  6  knts,  closest  approach  300  m) 


Figure  1  The  magnetic  dipole  field  of  a 
submarine  sailing  at  constant  speed  along  a 
straight  line,  as  seen  by  a  vector  magnetometer 
at  a  fixed  position,  is  always  and  in  all  its  vector 
components,  composed  of  linear  combinations 
of  three  simple  algebraic  functions  named 
** Anderson  functions The  coefficients  of  the 
linear  compositions  depend  on  the  dipole 
strength  and  orientation,  and  the  ship's  velocity, 
heading,  and  distance  from  the  magnetometer 
at  the  point  of  closest  approach. 


Montecristo  -  magnetic  field  vector  residuals 


OBM  04  ->  OBM  19  OBM  03  ->  OBM  19  OBM  04  ->  OBM  03 


Figure  2  Magnetic  field  vector 
residuals  for  three  pairs  of  OEMs, 
04  ->19,  03  ->  19  and  04  03 

(see  Map  2  for  their  positions). 
The  residuals  are  basically  the 
spatially  uncorrelated  components 
of  magnetic  field  variations 
observed  simultaneously  at  spaced 
sites.  The  relatively  large  excursion 
at  14:38  UTC  marks  the  passage 
of  the  workboat  Manning  across 
OEM  04.  See  the  text  for  further 
explanations. 


Figure  3  Mean  peak-to-peak  amplitudes  of 
natural  geomagnetic pc-type  variations  observed 
over  many  years  at  the  Geophysical  Institute 
Gottingen  (F.R.G.).  Such  variations  contribute  to 
the  always  present  geomagnetic  noise  and  set  a 
limit  to  the  submarine  detection  capability  if  no 
ambient  noise  suppression  were  applied. 
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Variation  of  the  ambient  magnetic  field 

observed  at  the  sea  bottom  by  OBM  04 


Figure  4a  Temporal  variation  of  the  three 
vector  components  of  the  ambient  magnetic 
north  on  a  magnetically  moderately 

disturbed  day  (upper  panel)  and  a 
magnetically  quiet  day  (lower  panel), 
recorded  with  OBM  04  (marked  *M04  **  in 
Map  3).  The  highlighted  wavelets  resemble 
vertical  in  shape,  amplitude  and  duration  the 
magnetic  field  distortion  generated  by  a 
08:40  08:45  08:50  08:55  mm  UJC  ship  (c.f.  Figure  1). 


date  95/09/16 


date  95/09/20 


Figure  4b  Variation  of  the  eastward 
component  of  the  ambient  magnetic  field 
recorded  simultaneously  with  three 
OEMs  at  spaced  sites.  The  time  intervals 
are  identical  to  those  shown  in  Figure  4a. 
OBM  18  and  OBM  19  were  located  2  km 
eastward  and  6  km  westward, 
respectively,  of  OBM  04  (see  also  Map 
3).  The  magnetic  field  difference  between 
OBM  04  and  the  mean  of  OBM  18  and 
19  (bottom  curve  in  each  of  the  two 
panels)  reveals  that  the  noted  wavelets 
were  uniform  over  several  kilometres 
distance  (upper panel)  and  localized  near 
OBM  04  (lower panel),  respectively.  The 
characteristics  of  the  mean  difference 
suggest  that  the  wavelets  were  of 
ionospheric  origin  in  the  first  case  (upper 
panel),  and  generated  by  a  passing  ship 
in  the  second  case  (lower  panel), 
respectively.  The  northward  and  vertical 
components  of  the  magnetic  variation 
were  in  principle  similar  to  the  eastward 
component  and  are  not  shown  here. 


Variation  of  the  ambient  magnetic  field 

observed  on  the  sea  bottom  at  spaced  sites 
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Figure  5  Diurnal  distribution  of  the  geomagnetic  field  variance  near  Deiva  Marina,  averaged  over 
consecutive  15-min  intervals.  Only  relative  amplitudes  are  displayed.  The  measurements  were  made 
continuously  over  16  days  with  a  land-based  magnetometer  operated  by  A.  Magunia  (Universitdt 
Frankfiirt/Main,  F.R.  G.)  at  about  500  m  distance  from  the  railway  line.  The  small  arrowheads  along 
the  top  and  bottom  mark  the  separations  we  chose  in  order  to  break  the  day  into  four  intervals  with 
different  levels  of  magnetic  field  variance  (c.f  Table  4). 


Ligurian  coast  1996 

Magnetic  field  component  perpendicular  to  coast 
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Figure  6  Magnetic  field  variations  (coast-perpendicular  component  only)  recorded  simultaneously 
by  the  OEMs  17, 18, 19  and  12,  which  were  aligned  along  a  profile  perpendicular  to  the  coast  and  to 
a  major  railway  line  (which  closely  follows  the  coast  line),  see  Map  1.  The  entire  time  interval  (1 600- 
1900  UTC,  equivalent  to  1 7:00-20:00  local  time)  is  strongly  affected  by  train  noise  which  significantly 
exceeds  the  ambient  magnetic  noise  from  natural  sources  even  at  8  km  distance  from  the  railway  line. 
Leakage  of  electric  train  currents  into  the  sea  is  the  most  likely  source  of  the  magnetic  noise. 


Time 


Bin  tfireshold 


Magnetic  field  variation 


O  Points  counted 


Figure  7  Schematic  of  the  binning 
method  used  to  create  the  distributions 
presented  in  Figures  8, 10  and  11.  Each 
time  the  magnetic  field  variation, 
difference,  or  residual  crosses  a  bin 
threshold  in  the  direction  away  from 
the  zero  line,  the  counter  is  incremented. 
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Figure  8  Occurrence  frequency  distributions  of the  magnetic field  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  C*events/hour**),  cf  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8a  refers  to  the  1997  data  from  the  magnetically  quiet  Montecristo  area. 
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Figure  8  Occurrence  frequency  distributions  of  the  magnetic  field  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  f 'events/hour'*),  c.f  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f.  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8b  refers  to  the  1997  data  from  the  magnetically  quiet  Montecristo  area. 
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Figure  8  Occurrence  frequency  distributions  of  the  magnetic  field  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  f 'events/hour  ),  cf.  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8c  refers  to  the  1996  data  from  the  highly  disturbed  southern  Ligurian  coastal  zone. 
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Figure  8  Occurrence  frequency  distributions  of  the  magnetic  field  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  (“events/hour'*),  c.f  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having four  diurnal  intervals,  cf  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8d  refers  to  the  1996  data  fivm  the  highly  disturbed  southern  Ligurian  coastal  zone. 
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Figure  8  Occurrence  frequency  distributions  of  the  magnetic  field  variations,  differences  and  residuals 
between  two  neighbouring  trUaxial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  events/hour**),  c.f  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f.  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8e  refers  to  the  1996  data  from  the  highly  disturbed  southern  Ligurian  coastal  zone. 
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Figure  8  Occurrence frequency  distributions  of the  magnetic  field  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  f'events/hour*f,  c.f  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8f  refers  to  the  1996  data  from  the  highly  disturbed  southern  Ligurian  coastal  zone. 
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Figure  8  Occurrence  frequency  distributions  of  the  magnetic  field  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  f  events/hour''),  c.f  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f.  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8g  refers  to  the  1996  data  from  the  highly  disturbed  southern  Ligurian  coastal  zone. 
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Figure  8  Occurrence  frequency  distributions  of  the  magnetic  freld  variations,  differences  and  residuals 
between  two  neighbouring  tri-axial  magnetometers.  The  OBM  identification  numbers  appear  on  the 
right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to  number  of  bin  threshold  crossings 
per  hour  f 'events/hour*'),  c.f  Figure  7.  The  bin  thresholds  increase  along  the  ordinate,  the  days 
(each  having  four  diurnal  intervals,  c.f.  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into 
two  subpanels,  the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  Tables  7  and  8). 
Figure  8h  refers  to  the  1996  data  from  the  highly  disturbed  southern  Ligurian  coastal  zone. 
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Magnetic  field  variation  during  train  passage 
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Figure  9  Variation  of  the  ambient  magnetic  field  (coast-parallel 
component),  measured  simultaneously  with  four  OEMs  located  along 
a  line  perpendicular  to  the  coast  (modules  17, 18, 19  and  12,  see  Map 
1).  The  highlighted  section  marks  the  passage  of  a  single  train  on  the 
coastal  railway  line.  Yellow  covers  the  interval  when  the  train  was  on 
the  rails  between  the  power  substations  ofLevanto  and  Sestri  Levante, 
green  the  passing  of  the  train  through  the  railway  station  of  Deiva 
Marina  (without  making  a  stop). 
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a) 


Montecristo:  OBM  17  vs.  reference  OBM  19  (part  I) 
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Figure  10  Cumulative  occurrence  frequency  distributions  of  magnetic  field  variations,  differences 
and  residuals  between  two  neighbouring  sites  (1997  data,  Montecristo  area).  Each  panel  displays 
the  results  from  a  specific  vector  component  and  a  specific  diurnal  interval,  averaged  over  the  length 
of  the  entire  sea  trial  The  term  ''events/hour  ”  has  the  same  meaning  as  in  Figure  8:  The  bar  diagram 
indicates  how  often  a  certain  bin  theshold,  the  size  of  which  increases  along  the  abscissa,  is  exceeded. 
The  comparison  between  OEMs  17  and  19  is  split  into  two  parts,  (I)  and  (U),  in  accordance  with 
Figure  8a.  10a  contains  the  data  from  the  first  part  of  the  sea  trial,  up  to  July  23,  and  10b  the  much 
noisier  data  from  the  second  part,  starting  July  23. 
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Montecristo:  OBM  17  vs.  reference  OBM  19  (part  II) 
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Figure  10  Cumulative  occurrence  frequency  distributions  of  magnetic  field  variations,  differences 
and  residuals  between  two  neighbouring  sites  (1997  data,  Montecristo  area).  Each  panel  displays 
the  results  from  a  specific  vector  component  and  a  specific  diurnal  interval,  averaged  over  the  length 
of  the  entire  sea  trial.  The  term  ‘'events/hour  ”  has  the  same  meaning  as  in  Figure  8:  The  bar  diagram 
indicates  how  often  a  certain  bin  theshold,  the  size  of  which  increases  along  the  abscissa,  is  exceeded. 
The  comparison  between  OEMs  17  and  19  is  split  into  two  parts,  (I)  and  (U),  in  accordance  with 
Figure  8a.  10a  contains  the  data  from  the  first  part  of  the  sea  trial,  up  to  July  23,  and  10b  the  much 
noisier  data  from  the  second  part,  starting  July  23. 
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Figure  10  Cumulative  occurrence  frequency  distributions  of  magnetic  field  variations,  differences 
and  residuals  between  two  neighbouring  sites  (1997  data,  Montecristo  area).  Each  panel  displays 
the  results  from  a  specific  vector  component  and  a  specific  diurnal  interval,  averaged  over  the  length 
of  the  entire  sea  trial  The  term  ''events/hour  **  has  the  same  meaning  as  in  Figure  8:  The  bar  diagram 
indicates  how  often  a  certain  bin  theshold,  the  size  of  which  increases  along  the  abscissa,  is  exceeded. 
The  comparison  between  OBMs  17  and  19  is  split  into  two  parts,  (I)  and  (II),  in  accordance  with 
Figure  8a.  10c  covers  the  entire  length  of  the  sea  trial 
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Figure  11  Like  Figure  10  but  for  the  1996  data  (southern  Ligurian  coast).  In  this  case,  the  bin  sizes 
were  normalized  to  a  hypothetical  sensor  spacing  of  1  km,  taking  into  account  the  non-normalizable 
dual-sensor  system  noise  (see  Annex  D  for  a  detailed  discussion).  Therefore  the  abscissa  does  not 
start  with  0  nT/km  but  with  a  certain  offset  which  varies  with  the  distance  between  the  sensors  in 
each  pair. 
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b) 


Ligurian  coast:  OBM  19  vs.  reference  OBM  18 
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Figure  11  Like  Figure  10  but  for  the  1996  data  (southern  Ligurian  coast).  In  this  case,  the  bin  sizes 
were  normalized  to  a  hypothetical  sensor  spacing  of  1  km,  taking  into  account  the  non-normalizable 
dual-sensor  system  noise  (see  Annex  D  for  a  detailed  discussion).  Therefore  the  abscissa  does  not 
start  with  0  nT/km  but  with  a  certain  offset  which  varies  with  the  distance  between  the  sensors  in 
each  pair. 
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Ligurian  coast:  OBM  12  vs.  reference  OBM  19 
X  component  Y  component  Z  component 
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Figure  11  Like  Figure  10  but  for  the  1996  data  (southern  Ligurian  coast).  In  this  case,  the  bin  sizes 
were  normalized  to  a  hypothetical  sensor  spacing  of  1  km,  taking  into  account  the  non-normalizable 
dual-sensor  system  noise  (see  Annex  D  for  a  detailed  discussion).  Therefore  the  abscissa  does  not 
start  with  0  nT/km  but  with  a  certain  offset  which  varies  with  the  distance  between  the  sensors  in 
each  pair. 
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Ligurian  coast:  OBM  06  vs.  reference  OBM  12 
X  component  Y  component  Z  component 
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Figure  11  Like  Fig^re  10  but  for  the  1996  data  (southern  Ligurian  coast).  In  this  case,  the  bin  sizes 
were  normalized  to  a  hypothetical  sensor  spacing  of  1  km,  taking  into  account  the  non-normalizable 
duaUsensor  system  noise  (see  Annex  D  for  a  detailed  discussion).  Therefore  the  abscissa  does  not 
start  with  0  nT/km  but  with  a  certain  offset  which  varies  with  the  distance  between  the  sensors  in 
each  pair. 
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e) 


Ligurian  coast:  OBM  04  vs.  reference  OBM  06 
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Figure  11  Like  Figure  10  but  for  the  1996  data  (southern  Ligurian  coast).  In  this  case,  the  bin  sizes 
were  normalized  to  a  hypothetical  sensor  spacing  of  1  km,  taking  into  account  the  non-normalizable 
dual-sensor  system  noise  (see  Annex  D  for  a  detailed  discussion).  Therefore  the  abscissa  does  not 
start  with  0  nT/km  but  with  a  certain  offset  which  varies  with  the  distance  between  the  sensors  in 
each  pair. 
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Mean  time  interval  between  magnetic  field 
fluctuations  with  residuals  >  0.4  nT/km 
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Figure  12  Mean  lifetime  of  magnetic  field  residuals  between  two  neighbouring  OEMs 
before  they  exceed  0.4  nT/km.  The  relevant  OEM  numbers  are  printed  along  the  bottom 
of  the  figure.  The  threshold  of  0.4  nT/km  includes  a  distance-independent  contribution  of 
0.25  nf  our  estimate  of the  maximum  dual-sensor  system  noise.  The  bar  diagram  suggests 
short  intervals  of  uninterrupted  low  residual  amplitudes  for  the  vertical  magnetic  field 
component  and  long  intervals  for  residuals  between  the  coast-perpendicular  component 
from  two  magnetometers  at  the  same  distance  from  the  coast. 
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Annex  A 

Transformation  of  a  vector  from  the 
magnetometer  coordinate  system  into  a 
geomagnetic  coordinate  system 


Let  us  assume  that  the  three  axes  of  a  vector  magnetometer  form  an  orthogonal, 
right-handed  coordinate  system  defined  by  the  unit  vectors  y^^^^and 

aligned  with  the  corresponding  sensor  axes.  It  is  our  objective  to  find  a  transfor¬ 
mation  of  the  magnetic  field  variation  vector,  B,  from  the  «(”*)-y("*)-i(”*)-system 
(i.e.  the  coordinate  system  in  which  the  components  of  B  were  measiured)  into  a 
geomagnetic  north,  east  and  vertical  (N-E-V)  coordinate  system. 

The  inclination  angles  of  the  and  y("*)-axes  against  the  horizontal  plane.  Ox 

and  6y,  are  assumed  to  have  been  measured.  They  are  signed  positive  when  the 
associated  axes  point  upward.  The  angle  between  Nadir  and  i('")-axis,  0^,  has  not 
been  measured  but  can  be  computed  from  9x  and  6y.  Specifically,  we  require  that 
the  following  constraints  apply  to  the  inclination  angles 


-IT  12  <  0x,0y 

<  7r/2 

0  <  Oz 

<  7r/2 

(Al) 

0x7^0  V 

^  0 

Fulfillment  of  these  conditions  implies  that  the  magnetometer  module  is  never  turned 
upside  down.  The  heading  angle.  A,  measured  from  the  projection  of  the  geomag¬ 
netic  field  vector,  F,  on  the  -plane,  to  the  magnetometer  x-sensor  (aligned 

with  the  x(’")-axis  by  definition),  signed  positive  clockwise  around  the  ;?('"^-axis,  is 
assumed  to  have  been  recorded  as  well.  If  both  inclination  angles.  Ox  and  were 
equal  to  zero,  the  transformation  problem  would  be  trivial  and  solved  by  a  simple 
rotation  of  the  -plane  around  the  i(’")-axis  by  the  angle  —A. 

The  angles  Ox,  Oy  and  A  are  occasionally  confused  with  the  Euler  angles  [16],  and 
the  latter,  in  turn,  with  the  attitude  angles,  pitch,  roll  and  yaw.  Euler  angles  are 
easy  to  handle  and  result  in  a  simple,  straightforward  transformation  scheme.  In 
reality,  however,  it  is  rarely  possible  to  measure  the  Euler  angles  directly.  We  will 
see  further  below  that  Euler  angles  appear  also  in  our  algorithm,  in  fact,  our  angles 
Oz,  V’  and  ip  correspond  to  the  Euler  angles  0,—^  and  — in  the  terminology 
used  in  [16].  However,  it  takes  some  trigonometric  and  algebraic  manipulation  to 
determine  them.  This  is  described  in  the  following. 
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We  define  an  auxiliary  coordinate  system  based  on  the  unit  vectors  ei,  @2  and  §3, 
in  such  a  way  that  63  is  pointing  vertically  down  (i.e.  §3  =  V),  ei  is  aligned  with 
the  projection  of  on  the  horizontal  plane,  and  §2  completes  the  orthogonal, 
right-handed  system.  Note  that  @2  lies  in  the  horizontal  plane  but  is  generally  not 
aligned  with  the  projection  of  on  the  horizontal  plane,  except  for  =  0  or 
6y  =  0.  The  vector  can  be  represented  in  the  ei-e2-e3-system  as 


/  COS  Ox 

\ 

^(m)  ^ 

0 

y  -sin^a; 

/ 

{A2) 


The  first  and  the  second  components  result  fi:om  the  definition  of  Ox,  the  third  fi:om 
the  requirement  for  unit  length,  and  the  sign  is  required  because  of  the  upward 
orientation  of  opposite  to  63,  when  Ox  is  signed  positive.  The  vector 


II 

?  ?  ? 

(A3) 

is  found  from  the  identities 

0 

II 

? 

? 

("»)  =  _  tan  Ox  sin  Oy 

(A4) 

|5r("»)|  =  1 

=  ±^l-  sin^  Oy  /  cos2  Ox 

In  the  last  equation,  only  the  positive  root  is  valid,  otherwise  the  module  would 
have  been  turned  upsided  down.  The  geometrical  constraints  on  the  attitude  ensure 
sin^  6y  <  cos^  Ox-  The  representation  of  in  the  ei-e2-e3-system  is  found  from 


a. 

±i‘m)  ^  ^(m) 

— 

tan  0a:  \/  cos^  Ox  —  sin^  Oy 

= 

sin  Oy!  cos  Ox 

(A5) 

= 

J  cos^  Ox  —  sin^  Oy 
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From  =  cos  6z  (by  definition)  we  obtain 


cos^  62  =  cos^  9x  —  sin^  9y  =  1  —  sin^  $x  —  sin^  Oy 

sin^  6z  =  sin^  6x  +  sin^  6y 


(A6) 


We  recollect  the  representations  of  x("*\  and  z^”*)  in  the  ei-e2-e3-system 


/  COS0a;  \ 

^  —  tan  0a;  sin  0y  ^ 

^  tan  Ox  cos  02 

5^(m)  ^ 

II 

? 

0 

COsOz!  cos  0a: 

sin0y/  COS0X 

^  -sm9x  ) 

-sin0y  y 

^  cos  02 

(A7) 


The  inverse  transform,  namely  the  representation  of  the  unit  vectors  ei,  @2  and  @3, 
in  the  £(’”)-y("*)-i(’”^-system,  can  be  found  after  a  little  algebra 


ei  = 


COS^I  \ 

—  tan^aisin^j,  1  @2  = 

tan  cos  02  / 


0 

COS0z/ cos  01 
sin0y/cos0i 


(-sin  0a:  \ 
-sin0y  j 
cos  02  J 


(A8) 


The  unit  vector  along  the  projection  of  z^"*^  onto  the  ei-e2-plane  reads 


1 

sin  02 


^  tan  0j;  cos  02  \ 
sin  0J,/ cos  0a;  j 

0  J 


(A9) 


and  the  unit  vector  obtained  by  rotating  Zp^”*)  clockwise  by  7r/2  around  the 

ea-axis,  reads 


^(m)  _ 


1 

sin  02 


—  sin0y/cos0a:  \ 
tan  0a;  cos  02  j 

0  J 


(AlO) 


Zp("*)  and  are  orthogonal  vectors,  both  confined  to  the  ei-e2-plane.  The 

angle,  ip,  between  and  w^"*^  measured  positive  clockwise  from  y(”*)  to 
is  determined  from  the  following  relations 
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cosip  =  yC*")  •  =  sin^x/sin^z 

sinV’-z^"*)  =  y(”*)  x  =>■  s'mip  =  s\n6y  j  sin^z 


(All) 


The  latter  equation  follows  from  the  fact  that  y^"*)  x  is  (i)  perpendicular  to 
which  confines  the  cross  product  to  the  -plane,  and  (ii)  perpendicular  to 

which  confines  the  cross  product  to  the  z(’^^-Zp(”*^-plane.  Both  planes  have 
the  vector  in  common. 

We  rotate  the  x("*)-y(”*)-plane  clockwise  around  the  z(”*)-axis  by  the  angle  ip  in 
order  to  construct  a  new,  right-handed  coordinate  system  based  on  the  orthonormal 
vectors,  and  z^"*^  where  xz^”*^.  In  the  new  system,  the  angle 

A  has  to  be  replaced  by  A*  =  A  + ip.  In  the  northern  hemisphere,  the  unit  vector  of 
the  quasi-static  geomagnetic  field,  F,  points  into  the  halfspace  below  the  horizontal 
plane.  F  can  be  expressed  through  the  northward  and  downward  pointing  unit 
vectors,  N  and  V,  with  the  help  of  the  geomagnetic  inclination  angle,  I: 

F  =  N  cos  /  -I-  V  sin  J  (A12) 


V  caji  be  expressed  through  a  combination  of  z^”*)  and 


V  =  cos^zz^"*^ -sin^zU^”*) 

4  (A13) 

F  =  N  cos  /  -I-  z^”*^  cos  6z  sin  J  —  sin0z  sin/ 

The  unit  vector  of  the  projection  of  F  upon  the  u^"‘^-ti)(”*)-plane  reads: 

fp^rr.)  =  COS  A*  -  sin  A*  (A14) 

Hence,  F  takes  the  form: 

F  =  a  z("^)  -t-  b  cos  A*  -  sin  A*)  (A15) 

with  yet  unknown  coefficients  a  and  b.  Obviously,  a  and  b  must  always  be  real,  and 
6  must  be  non-negative,  otherwise  eq.  (A15)  would  contradict  eq.  (A14).  We  obtain 
the  relation  between  a  and  b  from 
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1  =  F2  =  +  b'^icos^  A*  +  sin2  A*)  =  0^  +  62 

a  =  ±\/ 1  —  6^ 


(A16) 


In  order  to  simplify  the  following  calculations  we  limit  the  module  inclination,  0^, 
to  reasonable  angles.  We  conducted  our  measurements  in  the  northern  hemisphere, 
in  an  ztfea  where  /  «  60°,  and  we  require  0,  <  I.  The  latter  guarantees  a  >  0,  i.e. 
only  the  positive  root  is  valid.  In  fact,  in  all  our  sea  trials,  dz  was  typically  an  order 
of  magnitude  smaller  than  I.  Furtheron,  the  following  relations  hold 


F  •  =  b  cos  A* 

F.-^(”»)  =  —6  sin  A*  (A-17) 

F  •  z("»)  =  a 

Prom  eq.  (A13)  we  find 

F-ul”*)  =  N •  cos/ —  sin^xsin/ 

F-w("*)  =  N-w("*)cos/  (A18) 

F  •  =  N  •  z^”*)  cos  I  +  cos  dz  sin  / 


On  the  right-hand  side  of  eq.  (A17),  we  express  a  through  b  using  the  result  firom 
eq.  (A16).  We  then  combine  eq.  (A17)  with  eq.  (A18)  in  order  to  eliminate  the 
vector  products  involving  F.  Consider  now 


1  =  N2  =  (N  •  +  (N  •  +  (N  •  (A19) 


We  replace  the  three  squared  terms  on  the  right-hand  side  of  eq.  (A19)  with  the 
corresponding  terms  obtained  firom  eq.  (A18).  After  a  little  algebra  we  arrive  at  the 
quadratic  equation: 


62  +  2Ab-C^  =  0 


(A20) 


in  which  we  used  the  abbreviations 


A  = 


smIsm$zCOsA* 

sin2  02  cos^  A*  +  cos^  0^ 


C^  = 


cos2  02  —  sin2 1 

sin2  02  cos^  A*  +  cos^  02 


(A21) 
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If  we  require  +  /  <  7r/2,  a  condition  which  was  always  fulfilled  in  our  experiments, 
we  find  that  is  always  positive.  The  general  solution  reads 


b=-A±  +  (A22) 

Evaluation  of  the  relation 


2  ^2  _  Oz{cos^  A*  sin^  6z  +  cos^  6^  —  sin^  >  n 

^  (sin^  9z  cos^  A*  +  cos^  Oz)^  ~ 


(A23) 


guarantees  that  b  is  real.  As  b  must  be  non-negative,  only  the  solution  of  eq.  (A22) 
with  the  positively  signed  second  term  is  valid.  The  magnetic  variation  unit  vector, 
B,  and  the  geomagnetic  field  unit  vector,  F,  read  in  the  u(’")-tt;("*)-i('")-system 


I 

■S^(m)  / 

/  ^«(m)  \  /  6cosA*  \ 

=  -bsinA*  (A25) 

V  A,.,  /  V  / 


cosV’  sinV’  0 
—  sin^  cos‘ip  0 
0  0  1 


(A24) 


We  rotate  the  -system  by  —9z  around  the  «)("*)-axis,  in  order  to  con¬ 

struct  a  new  coordinate  system,  named  the  tt*-u)(”*)-.S*-system.  Note  that  z*=  €3, 
i.e.  the  -plane  is  the  horizontal  plane.  The  magnetic  variation  unit  vector, 

B,  reads  in  the  u*^("*)-i*-system: 


Bu>  N  /  cos  02  0  sin  02  \  /  B^(m)  \ 

=  0  1  0  (A26) 

Bz-  )  \  -sin 02  0  cos 02  j  \  B^(m)  j 


The  geomagnetic  field  unit  vector,  P,  is  formed  the  same  way,  using  the  same  trans¬ 
formation  matrix.  With  the  help  of  eq.  (A12),  the  projection  of  F  upon  the  u*- 
t&l"*) -plane,  Fp*  (which  is  not  a  unit  vector),  can  be  expressed  as 

Fp.  =  NCOS/  =  (F„.u*  -t-  4(„.)w(”*))cos/  (A27) 


The  angle  between  Fp*  and  the  fi*-axis,  signed  positive  when  measured  clockwise 
around  the  i*-axis  from  Fp.  to  u*,  is  determined  by 
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cos  v?  =  N  •  u*  = 


F  •  u*  _  Fu‘  _  cos 6zF^(m)  +smdzF^(m) 
cos  I  cos  I  cos  I 

_  b  cos  9z  cos  A*  +  v/ 1  -  6^  gin  6z 
cos  J 


(A28) 


—  sin^ 


=  N  •  = 


F„(m)  -6sinA* 


cos  I  cos  / 


cos  I 


Finally,  we  rotate  the  u*-w^'^'>-z* -system  counterclockwise  by  ip  (clockwise  by  —ip) 
around  the  i*-axis  (=  ea-axis)  and  arrive  at  the  N-E-V-system.  The  following 
transform  permits  us  to  express  the  unit  vector  of  the  magnetic  variation,  B,  in  the 
latter  system. 


Bn  ^ 

^  cosip 

’-simp 

M 

f  By* 

= 

simp 

COS(p 

0 

Bv  ) 

[  0 

0 

ly 

Bz- 

(A29) 


is  the  desired  representation  of  B  in  the  geomagnetic  north-east-vertical  {N-E-V)- 
system. 

Occasionally,  a  fmrther  rotation  of  the  coordinate  system  around  the  vertical  axis  is 
performed  in  order  to  align  the  two  horizontal  axes  with  directions  imposed  by  the 
environment,  e.g.  the  topography  or  geology.  For  measurement  sites  close  to  the 
coast,  the  coast-parallel  and  coast-perpendicular  directions  are  an  obvious  choice. 
If  the  angle  from  magnetic  north  to  the  orientation  of  the  coast  line  is  denoted  by  7 
(measured  positive  clockwise  from  north),  a  further  matrix  transformation  needs  to 
be  applied  to  the  magnetic  field  vector. 


i 

cos  7 

sin  7 

(  Bn 

r 

—  sin  7 

cos  7 

0 

Be 

\Bv  )  \ 

.  0 

0 

1 J 

\Bv 

(A30) 


is  the  desired  representation  of  B  in  a  system  where  the  three  axes  are  oriented  in 
the  directions,  coast-parallel,  coast-perpendicular  and  vertical. 
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In  summary,  we  have  to  proceed  with  the  following  sequence  of  computations  in 
order  to  achieve  our  goal  and  arrive  at  a  representation  of  the  magnetic  variation 
vector,  B,  in  the  geomagnetic  north-east-vertical  system 


cosdz 

sinOz 
cosip 
sinip 
sin  A* 
cos  A* 

D 
A 

b 

COS(p 
SilKp 

After  these  computation  steps  have  been  completed  (only  once  for  every  magne¬ 
tometer  and  every  site,  i.e.  for  a  fixed  set  of  angles.  Ox,  Ox,  A)  we  have  to  transform 
each  B-vector  in  the  time  domain  for  each  time  datum  according  to  the  scheme 


(B||(t)  \  /  COS7  sinq  0  \  /  cosip  —sintp  0  \ 

B^it)  I  =  —  sin7  COS7  0  I  •  I  sin^  cos^?  0  I 

Bv{t)  )  \  0  0  I  )  \  0  0  1  ) 

cosOz  0  sin0^  \  /  cosip  sinip  0  \  \ 

0  1  0  I  •  —sinip  cosip  0  •  j  I 

-sin 02  0  cos 02  /  \  0  0  1  /  \  B^(m){t)  j 


=  y  sin^  Ox  +  sin^  Oy 
=  sin  0i/ sin  02 
=  sinOy!  sin  02 
=  sin  A  cos  Ip  -t-  cos  A  sin^ 

=  cos  A  cos  V’  —  sin  A  sin  V’ 

=  sin^  02  cos^  A*  +  cos^  Oz 
=  (sin  I  sin  Oz  cos  A*)  /  D 
=  (cos^  02  —  sin^  I)  /  D 

=  -a  +  Va^  +  c^ 

=  (6  cos  02  cos  A*  -h  V I  —  b^  sin  Oz)  /  cos  I 
=  6sinA*/  cos  I 
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Annex  B 

Estimating  the  coefficients  of  a 
linear  trivariate  model 

We  assume  that  a  discrete  time  series,  u{t),  can  be  expressed  through  a  linear 
combination  of  three  independent  discrete  time  series,  x{t),  y{t)  and  z{t),  plus  an 
uncorrelated  term,  5u{t): 

u{t)  =  ax{t)  + /3y{t) +'yz{t)  +  6u{i)  (Bl) 

The  residual  sum  of  squares  reads 

—  2  (ay^  ux  +  I^Y^uy  + ')Y^uz)  (B2) 

+  2  +  )97  XI 

Taking  partial  derivatives  yields  the  system 

^  -T^^y+^12^y+'yJ2y^] 

In  order  to  find  that  parameter  set,  a,  /3,  7,  which  minimizes  Q^,  we  solve  the  system 
of  three  normal  equations 
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dQ^ 

da 

df3 

&y 


0 

0 

0 


(B4) 


which  can  be  rewritten  as  a  system  of  three  linear  inhomogeneous  equations  with 
the  three  xmknown  variables,  a,  /3,  r- 


a^x"^  +  P^xy  +  'Y^zx  =  ^ux 

aY^xy  + +  'r'^yz  =  (B5) 

a  za:  +  /3  ^  y^r  +  7  ^ 

The  system  is  solved  using  standard  techniques  of  numerical  mathematics. 

The  same  procedure  is  applied  to  each  of  the  three  components  of  the  magnetic 
field  variation  vector,  where  u{t)  stands  for  any  of  the  test  magnetometer  vector 
components,  and  a:(<),  y{t)  and  z{t)  for  the  vector  components  of  the  reference 
magnetometer. 
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Annex  C 

Inverting  test  and  reference  sensor: 
The  effect  on  the  residual  sum  of  squares 


Let  us  assume  that  two  different  sensors,  with  identification  numbers  ”1”  and  ”2”, 
record  simultaneously  iV  samples  each  of  two  discrete  scalar  time  series,  vi{t)  and 
V2{t),  consisting  of  different  multiples  of  a  common,  time- varying  stochastic  signal, 
s{i),  and  Gaussian  white  noise,  ni(t)  and  n2(<),  respectively: 


vi(<)  =  Cl  s(t) -fni(t) 
V2(t)  =  C2  s{t)  +  n2{t) 


(Cl) 


Here,  ci  and  C2  are  real  constants  (scale  factors).  The  noise  may  stem  from  the 
system  and  from  external  noise  sources  which  are  uncorrelated  with  the  signal  and 
uncorrelated  between  b 
values  {i,j  e  {1,2}): 


sensors. 

Signal  and  noise 

E?{s} 

=  0 

E{s^} 

= 

E{ni} 

=  0 

EM 

= 

£?{sni} 

=  0 

E{ninj} 

=  0  for  z  ^  j 

(C2) 


The  conditions  listed  in  (C2)  state  that  the  noise  variance  is  assumed  to  be  identical 
at  both  sites.  Each  of  the  time  series  shall  be  expressed  through  the  other  in  a  linear 
univariate  representation  {i,j  G  {1,2}  and  i  ^  j): 


Vj(t)  =  aiVi{t)  +  5j{t) 


(C3) 


where  Sj  represents  that  part  of  vj  which  is  not  linearly  correlated  with  Uj.  The 
residual  sum  of  squares,  reads 
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Q/  =  =  Y^{vj-otiVif  (C4) 

In  order  to  find  the  minimum  value  of  Qj2  (and  thus  the  optimal  regression  coeffi¬ 
cient,  ai)  we  take  its  derivative  with  respect  to 


dai 


and  solve  the  corresponding  normal  equation 


dai 


=  0 


CX<i 


We  insert  into  eq.  (C4)  and  obtain 


Qj 


2 


E«i2 


Expansion  of  the  terms  on  the  right  hand  side  yields 


(C5) 


(C6) 


(C7) 


ViVj  =  CiCj  +  Ci  53  srij  +  Cj  53  + 13 

53  =  Cj^53^^+2cj53^”»+13”»^ 

13  V  =  c/53«^+2cj13«”j  +  I3”/ 

Taking  into  consideration  the  conditions  stated  in  eq.  (C2),  we  find  the  asymptotic 
limits  (i.e.  the  limits  for  N  -4  oo)  of  the  normalized  versions  of  the  three  terms  on 
the  left-hand  side  of  eq.  (C8) 

+  (C9) 
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SO  that  the  expectation  value  of  the  univariate  regression  coefficient,  Oj,  reads 

Rewriting  eq.  (C7)  with  the  normalized  terms  from  eq.  (C9)  yields  an  asymptotic 
limit  for  the  residual  sum  of  squares 


lim 

N-^CX> 


+  CTn^){Cj^(rs^  +  On^)  -  V/ 
a'^Os^  + 


(Cll) 


We  apply  the  same  procedure  to  and  find  the  asymptotic  limit  of  the  ratio 
between  and 


lim 

N-^oo 


Cj^(7s^  +  O-n^ 
Ci^as‘^  +  an^ 


(C12) 


eq.  (C12)  lends  support  to  the  following  interpretation:  If,  for  instance,  the  time 
series  V2  has  a  larger  variance  than  ui,  which  is  equivalent  to  requiring  that  the 
absolute  value  of  C2  exceeds  that  of  ci  (see  eq.  C9)),  the  ratio  of  the  residual  sums 
of  squares,  is  asymptotically  greater  than  1,  in  other  words,  exceeds 

Q\^.  That  means,  expressing  V2  linearly  through  vy  in  a  univariate  way,  c.f.  eq.  (C3), 
results  in  an  asymptotically  larger  error  than  expressing  wi  linearly  through  W2. 
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Annex  D 

Magnetometer  system  noise  and  its 
effect  on  the  residuals 


The  distribution  of  the  magnetic  field  residuals  between  neighbouring  magneto¬ 
meters  operated  in  a  magnetically  quiet  area  provides  us  with  an  upper  limit  of 
the  system  noise  over  the  bandwidth  considered  in  this  report,  1.5-80  mHz.  Based 
on  the  results  displayed  in  Fig.  10a,  we  assume  a  maximum  of  0.25  nT  for  the  sys¬ 
tem  noise,  et-r,  between  two  properly  functioning  magnetometers  denoted  by  the 
subscripts  T  and  R. 

If  all  magnetometers  axe  practically  identical  as  far  as  the  sensor  system  noise  vari¬ 
ance  is  concerned,  the  single  sensor  self  noise  is  found  according  to: 


4- 


R  —  4+4 


=  £R 


£T-R 


0.25 

s/2 


0.18  nT 


(Dl) 


We  conclude  that  the  system  noise  of  each  sensor  component  of  each  magnetometer 
is  less  than  0.2  nT  in  the  bandwidth  of  interest. 

Let  us  now  consider  a  magnetic  field  residual,  Aobs,  computed  firom  observations 
from  the  same  vector  component  of  two  magnetometers,  subscripted  T  and  R,  which 
are  separated  by  a  distance,  D.  The  occurrence  frequency  of  those  residuals  which 
exceed  a  certain  threshold  decreases  rapidly  with  increasing  threshold.  We  thus 
argue  that  the  vast  majority  of  observed  excesses  of  an  arbitrary  threshold  occurs 
when  large  system  noise  and  true  magnetic  field  residual  add  up. 


A  /»  1  „  _ Atrue  ■Ao6s  ^T-R 

Aobs  =  Atrue  +  £T-R  ^ -  \^2) 

If  we  want  to  know  what  magnetic  field  residual  would  have  been  observed  if  the 
magnetometer  spacing  were  normalized  to  a  unit  distance  and  the  true  residual  had 
a  constant  spatial  gradient  (i.e.  the  residual  were  proportional  to  the  distance),  we 
need  to  find  the  true  normalized  amplitude,  supplemented  by  a  contribution  firom 
the  sensor  system  noise  which  is  independent  of  the  sensor  spacing.  In  doing  so  we 
arrive  at  the  following: 
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^  +  £T-«  =  ^  +  £T-H  (l  -  i)  m 

According  to  this  formula,  the  observed  magnetic  field  residual,  Agba:  would  have  to 
exceed  a  threshold.  Lobs,  given  by 


Lobs  —  D  Lunit  ~  £t-r  (D  —  1)  (D4) 


in  order  to  find  a  imit-length  spaced,  but  noise  contaminated,  residual  (i.e.  the  left- 
hand  side  of  eq.  (D3)),  to  exceed  a  threshold,  Lunit-  The  numbers  for  Lobs  given  in 
Table  5  of  our  report  were  computed  for  a  threshold,  Lunit  —  0-4  nT,  using  eq.  (D4). 
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Map  1  The  southern  Ligurian  coastal  zone  with  magnetometer  sites  marked  by 
black  dots.  The  red  dot  in  the  inlay  map  of  Italy  indicates  the  geographic  location  of 
the  measurement  area.  Bathymetry  in  metres. 

Map  2  Montecristo  island  (Tyrrhenian  Sea)  and  its  easterly  water  space,  with 
magnetometer  sites  marked  by  black  dots.  The  red  dot  in  the  inlay  map  of  Italy 
indicates  the  geographic  location  of  the  measurement  area.  Bathymetry  in  metres. 

Map  3  The  Formiche  di  Grosseto  area  (Tyrrhenian  Sea)  with  part  of  Tuscany 
visible  in  the  upper  right  comer.  Magnetometer  sites  are  marked  by  black  dots.  The 
red  dot  in  the  inlay  map  of  Italy  indicates  the  geographic  location  of  the  measurement 
area.  Bathymetry  in  metres. 

Figure  1  The  magnetic  dipole  field  of  a  submarine  sailing  at  constant  speed  along 
a  straight  line,  as  it  is  seen  by  a  vector  magnetometer  at  a  fixed  location,  is  always 
and  in  all  its  vector  components,  composed  of  linear  combinations  of  three  sim¬ 
ple  algebraic  functions  named  ” Anderson  functions”.  The  coefficients  of  the  linear 
compositions  depend  on  the  dipole  strength  and  orientation,  and  the  ship^s  velocity, 
heading,  and  distance  from  the  magnetometer  at  the  point  of  closest  approach. 

Figure  2  Magnetic  field  vector  residuals  for  three  pairs  of  OBMs,  04  19,  03  19  and 
04  03  (see  Map  2  for  their  locations).  The  residuals  are  basically  the  spatially  un¬ 
correlated  components  of  magnetic  field  variations  observed  simultaneously  at  spaced 
sites.  The  relatively  large  excursion  at  14:38  UTC  marks  the  passage  of  the  workboat 
Manning  across  OBM  04-  See  the  text  for  further  explanations. 

Figure  3  Mean  peak-to-peak  amplitudes  of  natural  geomagnetic  pc-type  variations 
observed  over  many  years  at  the  Geophysical  Institute  Gottingen  (F,R.G.),  Such 
variations  contribute  to  the  always  present  geomagnetic  noise  and  set  a  limit  to  the 
submarine  detection  capability  if  no  ambient  noise  suppression  is  applied. 

Figure  4a  Temporal  variation  of  the  three  vector  components  of  the  ambient  mag¬ 
netic  field  on  a  magnetically  moderately  disturbed  day  (upper  panel)  and  a  magnet¬ 
ically  quiet  day  (lower  panel),  recorded  with  OBM  04  (marked  ”M04”  in  Map  3), 
The  highlighted  wavelets  resemble  in  shape,  amplitude  and  duration  the  magnetic 
field  distortion  generated  by  a  passing  ship  (c.f.  Figure  1), 

Figure  4b  Variation  of  the  eastward  component  of  the  ambient  magnetic  field 
recorded  simultaneously  with  three  OBMs  at  spaced  sites.  The  time  intervals  are 
identical  to  those  shown  in  Figure  OBM  18  qnd  OBM  19  were  located  2  km 
eastward  and  6  km  westward,  respectively,  of  OBM  04  (see  also  Map  3).  The  mag¬ 
netic  field  difference  between  OBM  04  and  the  mean  of  OBM  18  and  19  (bottom 
curve  in  each  of  the  two  panels)  reveals  that  the  noted  wavelets  were  uniform  over 
several  kilometres  distance  (upper  panel)  and  localized  near  OBM  04  (lower  panel), 
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respectively.  The  characteristics  of  the  mean  difference  suggest  that  the  wavelets  were 
of  ionospheric  origin  in  the  first  case  (upper  panel),  and  generated  by  a  passing  ship 
in  the  second  case  (lower  panel),  respectively.  The  northward  and  vertical  compo¬ 
nents  of  the  magnetic  variation  were  in  principle  similar  to  the  eastward  component 
and  are  not  shown  here. 

Figure  5  Diurnal  distribution  of  the  geomagnetic  field  variance  near  Deiva  Marina, 
averaged  over  consecutive  15-min  intervals.  Only  relative  amplitudes  are  displayed. 
The  measurements  were  made  continuously  over  16  days  with  a  land-based  magne¬ 
tometer  operated  by  A.  Magunia  (Universitdt  Frankfurt/Main,  F.R.G.)  at  about  500 
m  distance  from  the  railway  line.  The  small  arrowheads  along  the  top  and  bottom 
margins  mark  the  separations  we  chose  in  order  to  break  the  day  into  four  intervals 
with  different  levels  of  magnetic  field  variance  (c.f.  Table  4)- 

Figure  6  Magnetic  field  variations  (coast-perpendicular  component  only)  recorded 
simultaneously  by  the  OEMs  17,  18,  19  and  12  which  were  aligned  along  a  profile 
perpendicular  to  the  coast  and  to  a  major  railway  line  (which  closely  follows  the 
coast  line),  see  Map  1.  The  entire  time  interval  (1600-1900  UTC,  equivalent  to 
17:00-20:00  local  time)  is  strongly  affected  by  train  noise  which  significantly  exceeds 
the  ambient  magnetic  noise  form  natural  sources  even  at  8  km  distance  from  the 
railway  line.  Leakage  of  the  electric  train  return  currents  into  the  sea  is  the  most 
likely  source  of  the  magnetic  noise. 

Figure  7  Schematic  of  the  binning  method  used  to  create  the  distributions  presented 
in  Figures  8,  10  and  11.  Each  time  the  magnetic  field  variation,  difference,  or 
residual  crosses  a  bin  threshold  in  a  direction  away  from  the  zero  line,  the  counter 
is  incremented. 

Figure  8  Occurrence  frequency  distributions  of  the  magnetic  field  variations,  dif¬ 
ferences  and  residuals  between  two  neighbouring  tri-axial  magnetometers.  Each  pair 
of  neighbouring  OEMs  occupies  one  printed  page.  The  OEM  identification  numbers 
appear  on  the  right  hand  side.  The  logarithmic  scale  at  the  bottom  relates  colour  to 
number  of  bin  threshold  crossings  per  hour  (” events/hour”),  c.f.  Figure  7.  The  bin 
thresholds  increase  along  the  ordinate,  the  days  (each  having  four  diurnal  intervals, 
c.f.  Table  4)  along  the  abscissa.  Each  residual  panel  is  split  into  two  subpanels, 
the  upper  one  containing  all  data,  the  lower  one  the  data  cleaned  for  identified  ship 
signatures,  acoustic  interrogation  noise,  and  certain  types  of  sensor  noise  (see  also 
Tables  7  and  8).  Figures  8a,b  refer  to  the  1997  data  from  the  magnetically  quiet 
Montecristo  area.  Figures  8c-h  to  the  1996  data  from  the  highly  disturbed  southern 
Ligurian  coastal  zone. 

Figure  9  Variation  of  the  ambient  magnetic  field  (coast-parallel  component),  mea¬ 
sured  simultaneously  with  four  OEMs  located  along  a  line  perpendicular  to  the  coast 
(modules  17,  18,  19  and  12,  see  Map  1).  The  highlighted  section  marks  the  passage 
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of  a  single  train  on  the  coastal  railway  line.  Yellow  covers  the  interval  when  the 
train  was  on  the  rails  between  the  power  substations  of  Levanto  and  Sestri  Levante, 
green  the  passing  of  the  train  through  the  railway  station  of  Deiva  Marina  (without 
making  a  stop). 

Figure  10  Cumulative  occurrence  frequency  distributions  of  magnetic  field  varia¬ 
tions,  differences  and  residuals  between  two  neighbouring  sites  (1997  data,  Monte- 
cristo  area).  Each  panel  displays  the  results  from  a  specific  vector  component  and  a 
specific  diurnal  interval,  averaged  over  the  length  of  the  entire  sea  trial.  The  term 
’’events/hour”  has  the  same  meaning  as  in  Figure  8:  The  bar  diagram  indicates 
how  often  a  certain  bin  theshold,  the  size  of  which  increases  along  the  abscissa,  is 
exceeded.  The  comparison  between  OEMs  17  and  19  is  split  into  two  parts,  (I)  and 
(II),  in  accordance  with  Figure  8a.  10a  contains  the  data  from  the  first  part  of  the 
sea  trial,  up  to  July  23,  and  10b  the  much  noisier  data  from  the  second  part,  starting 
July  23.  10c,  the  comparison  between  OEMs  04  and  19,  covers  the  entire  length  of 
the  sea  trial. 

Figure  11  Like  Figure  10  but  for  the  1996  data  (southern  Ligurian  coast).  For 
this  figure,  the  bin  sizes  were  normalized  to  a  hypothetical  sensor  spacing  of  1  km, 
taking  into  account  the  non-normalizable  dual-sensor  system  noise  (see  Annex  D  for 
a  detailed  discussion).  Therefore  the  abscissa  does  not  start  with  0  nT/km  but  with 
a  certain  offset  which  varies  with  the  distance  between  the  sensors  in  each  pair. 

Figure  12  Mean  lifetime  of  magnetic  field  residuals  between  two  neighbouring 
OEMs  before  they  exceed  O.4  nT/km.  The  relevant  OEM  numbers  are  printed  along 
the  bottom  of  the  figure.  The  threshold  of  0. 4  nT/km  includes  a  distance-independent 
contribution  of  0.25  nT,  which  is  our  estimate  for  the  dual-sensor  system  noise.  The 
bar  diagram  suggests  short  intervals  of  uninterrupted  low  residual  amplitudes  for  the 
vertical  magnetic  field  component  and  long  intervals  for  residuals  between  the  coast- 
perpendicular  component  from  two  magnetometers  having  the  same  distance  from 
coast. 
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Abstract 


Temporal  and  spatial  variations  of  the  geomagnetic  field  were  recorded  in  different  geographic  areas,  using  arrays 
of  seven  tri-axial  magnetometers  which  resided  on  the  sea  bottom  in  50-150  m  depth  for  periods  up  to  16  days. 
The  different  geographic  areas  are  characterized  by  different  levels  of  the  mean  ambient  magnetic  noise.  We 
discuss  briefly  data  from  a  magnetically  quiet  area  and  then  focus  on  a  coastal  zone  highly  disturbed  by 
anthropogenic  magnetic  fields,  in  particular  by  noise  from  a  busy  coastal  electrified  railway  line.  In  the  latter  area, 
the  spatial  uniformity  of  the  ambient  magnetic  field  is  rather  poor  and  the  correlation  between  adjacent  sites  lower 
than  in  the  former  one. 

Simultaneous  magnetic  field  observations  from  pairs  of  neighbouring  magnetometers  are  analysed  and  compared 
in  detail.  We  compute  first  vector  differences  of  the  magnetic  variations  measured  at  neighbouring  sites,  and 
subsequently  vector  residuals.  The  residuals  are  those  contributions  to  the  magnetic  field  variations  which  are  not 
correlated  betwedtt  adjacent  site  and  which  can  not  be  represented  through  a  linear  trivariate  model. 

In  the  magnetically  quiet  area,  we  find  that  the  magnetic  field  residuals  are  of  the  order  of  the  system  no«^^^^ 
magnetometers  with  about  1  km  spacing  and  slightly  higher  fpr  tho^  with  12  km  sp^hg^The  a^|l^4^;^j^l^iS  • 
to  the  application  of  a  Remote  Reference  Technique,  i.e.  a  technique  in  which  the  ambient  is 

compensated  for  by  using  measurements  from  a  remotely  operated  ;rcfeience  j^hsiof^\^^^^^  means  that  the 

reference  sensor  is  sufficiently  far  away  so  that  it  is  not  influenced  by  the  itti^gnelic  field  of  a  ship  navigating  in 
the  vicinity  of  the  surveillance  magnetom^^. 

Although  the  magnetic  field  variations  in  thd  hij^y  disturbed  coastal  zone  are  several  times  larger  than  typical 
magnetic  field  variations  of  natural  origin,  and  also  spatially  npn-uniform,  we  find  good  linear  correlation 
between  adjacent  sites.  This  is  demonstrated  by  the  often  rather  sni^l  magnetic  field  residuals.  However,  in  order 
to  achieve  small  residuals,  the  spacing  between  the  magnetometer  must  be  an  order  of  magnitude  sfhaller  than  in 
the  quiet  zone.  This  may  require  the  use  of  a  Local  Referent^  Technique  for  compensation  bf  the  ambient 
magnetic  field.  In  a  Local  Reference  Technique,  the  reference  sensor  must  be  so  close  to  the  surveillance  sensors 
that  it,  too,  is  influenced  by  the  magnetic  field  of  the  submarine  navigating  in  the  vicinity  of  the  surveillance 
sensors. 

In  the  coastal  zone,  we  find  that  the  residuals  of  the  vertical  component  of  the  magnetic  field  variations  remain 
rather  large.  The  coast-perpendicular  component  of  magnetometers  at  equal  distance  from  the  coast  yields  the 
smallest  residuals,  i.e.  they  yield  the  best  compensatipn  of  magnetic  field  variations  when  using  a  reference 
technique.  The  intensity  distribution  of  magnetic  field  variations  observed  at  the  various  sites  strongly  suggests 
that  a  significant  part  ol  the  electric  leakage  currents  from  the  railway  system  flows  out  into  the  sea,  up  to  several 
kilometres  off  the  shore. 

In  summary,  we  find  that  magnetic  field  variations  in  the  southern  Ligurian  coastal  zone  are  dominated  by 
magnetic  railway  noise  which  results  in  an  anisotropic  distribution  of  the  spatial  correlation  between  spaced  sites. 
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